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NYO-7689 

ON  THE  PARTIAL  DIFFERENCE  EQUATIONS  OF 
MATHEMATICAL  PHYSICS 

R.  Courant,  K.  Frledrlchs,  and  H.  Lewy 

Introduction 
Problems  Involving  the  classical  linear  partial  dif- 
ferential equations  of  mathematical  physics  can  be  reduced 
to  algebraic  ones  of  a  very  much  simpler  structure  by 
replacing  the  differentials  by  difference  quotients  on  some 
(say  rectilinear)  mesh.   This  paper  will  undertake  an  ele- 
mentary discussion  of  these  algebraic  problems.  In  partic- 
ular of  the  behavior  of  the  solution  as  the  mesh  width 
tends  to  zero.   For  present  purposes  we  limit  ourselves 
mainly  to  simple  but  typical  cases,  and  treat  them  In 
such  a  way  that  the  applicability  of  the  method  to  more 
general  difference  equations  and  to  those  with  arbitrarily 
many  independent  variables  Is  made  clear. 

Corresponding  to  the  correctly  posed  problems  for 
partial  differential  equations  we  will  treat  boundary  value 
and  eigenvalue  problems  for  elliptic  difference  equations, 
and  Initial  value  problems  for  the  hyperbolic  or  parabolic 
cases.   We  will  show  by  typical  examples  that  the  passage 
to  the  limit  Is  Indeed  possible,  I.e.  that  the  solution  of 
the  difference  equation  converges  to  the  solution  of  the 
corresponding  differential  equation;  In  fact  we  will  find 
that  for  elliptic  equation  In  general  a  difference  quotient 
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NYO-7689 
of  arbitrarily  high  order  tends  to  the  corresponding  deriva- 
tive.  Nowhere  do  we  assume  the  existence  of  the  solution 
to  the  differential  equation  problem  —  on  the  contrary, 
we  obtain  a  simple  existence  proof  by  using  the  limiting 
process  .   For  the  case  of  elliptic  equations  convergence 
is  obtained  independently  of  the  choice  of  mesh,  but  we 
will  find  that  for  the  case  of  the  initial  value  problem 
for  hyperbolic  equations,  convergence  is  obtained  only  if 
the  ratio  of  the  mesh  widths  in  different  directions  satis- 
fies certain  inequalities  which  in  turn  depend  on  the  position 
of  the  characteristics  relative  to  the  mesh. 

We  take  as  a  typical  case  the  boundary  value  problem 
of  potential  theory.   Its  solution  and  Its  relation  to  the 

solution  of  the  corresponding  difference  equation  has  been 

2 
extensively  treated  during  the  past  few  years.    However 

in  contrast  to  the  present  paper,  the  previous  work  has 


Our  method  of  proof  may  be  extended  without  difficulty  to 
cover  boundary  value  and  eigenvalue  problems  for  arisltrary 
linear  elliptic  differential  equations  and  initial  value 
problems  for  arbitrary  linear  hyperbolic  differential 
equations. 

J.  leRoux,  "Sur  le  problem  de  Dlrlchlet",  Journ.  de  mathem. 
pur.  et  appl.  (6)  10  (191^) »  P-  189.  R.O.D.  Richardson,  "A 
new  method  In  boundary  problems  for  differential  equations". 
Trans,  of  the  Am.  Math.  Soc.  I8  (1917) >  P-  ^89  ff.  H.B. 
Philips  and  N.  Wiener.  "Nets  and  the  Dlrlchlet  Problem", 
Publ.  of  M. I.T.  (1925). 

Unfortunately  these  papers  were  not  known  by  the  first  of 
the  three  authors  when  he  prepared  his  note  "On  the  theory 

-  2  - 


NYO-7689 
Involved  the  use  of  quite  special  characteristics  of  the 
potential  equation  so  that  the  applicability  of  the  method 
used  there  to  other  problems  has  not  been  Immediately 
evident. 

In  addition  to  the  main  part  of  the  paper,  we  append 
an  elementary  algebraic  discussion  of  the  connection  of 
the  boundary  value  problem  of  elliptic  equations  with  the 
random  walk  problem  arising  in  statistics. 


of  partial  difference  equations",  GiJtt.  Nachr.  23,  X,  1925, 
on  which  the  present  work  depends. 

See  also  L.  Lusternik,  "On  an  application  of  the  direct 
method  in  variation  calculus",  Recuell  dela  Societe' 
Ma  them,  de  Moscou,  I926.   G.  Bouligand,  "Sur  le  probleme 
de  Dlrichlet,"  Ann.  de  la  soc.  polon.  de  math^ra.  4^, 
Cracow  1926. 

On  the  meaning  of  the  difference  expressions  and  on  further 
applications  of  them,  see  R.  Courant,  "Uber  dlrekte  Methoden 
In  der  Variationsrechnung",  Math.  Annalen  97>  P-  711 >  and 
the  references  given  therein. 
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I.   The  Elliptic  Case 

Section  1.   Preliminary  Remarks 

1.   Definitions 

Consider  a  rectangular  array  of  points  In  the 
(x,y) -plane,  such  that  for  mesh  width  h  >  0  the 
points  of  the  lattice  are  given  by 

X  =  nh  ") 

>    m,n  =  0,+l,  +2,* • • . 
y  «  mh   3 

Let  Q  be  a  region  of  the  plane  bounded  by  a 
continuous  closed  curve  which  has  no  double  points.   Then 
the  corresponding  mesh  region,  Q^  —  which  Is  uniquely 
determined  for  sufficiently  small  mesh  width  —  consists 
of  all  those  mesh  points  lying  in  G  which  can  be 
connected  to  any  other  given  point  in  G  by  a  connected 
chain  of  mesh  points.  By   a  connected  chain  of  mesh  points 
we  mean  a  sequence  of  points  such  that  each  point  follows 
in  the  sequence  one  of  its  four  neighboring  points.   We 
denote  as  a  boundary  point  of  G.   a  point  whose  four 
neighboring  points  do  not  all  belong  to  G.  .   All  other 
points  of  G.   we  call  interior  points. 

We  shall  consider  functions  u,v,'''  of  position  on 
the  grid,  i.e.  functions  which  are  defined  only  for  grid 
points,  but  we  shall  denote  them  as  u(x,y) ,v(x,y) , • • • . 
For  their  forward  and  backward  difference  quotients  we 
employ  the  following  notation, 
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u^  =  ^[u(x+h,y)-u(x,y)],   u^  =  •p^[u(x,y+h)-u(x,y)  ]   . 

u-  =  ^[u(x,y)-u(x-h,y)],   uy  =  j^[u(x,y)-u(x,y-h)  ]   . 

In  the  same  way  the  difference  quotients  of  higher 
order  are  formed,  e.g. 

(^x^-  "  ^xx  "  ^x  "^   ^  [u(x+h,y)-2u(x,y)+u(x-h,y)],  etc. 
X  h 

2.   Difference  Expressions  and  Green's  Function 
In  order  to  study  linear  difference  expressions  of 
second  order,  we  form  (using  as  a  model  the  theory  of 
partial  differential  equations) ,  a  bilinear  expression 
from  the  forward  difference  quotients  of  two  functions, 
u  and  V, 

B(u,v)  =  auj^v^  +  ^^x^y  ■*■  ^Vx  "*"  ^Vy  "*"  °'"x^  "*■ 

3u  V  +  yuv  +  5uv  +  guv 
y       X      y 

where 

a  =  a(x,y),"-,  a  =  a(x,y),---,  g  =  g(x,y) 

are  functions  on  the  mesh. 

Prom  the  bilinear  expression  of  first  order  we  derive 
a  difference  expression  of  second  order  in  the  following 
way:  we  form  the  sura 

h^l     Y.   B(u,v) 
^h 
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over  all  points  of  a  region  G.   In  the  mesh,  where  In 
B(u,v)  the  difference  quotients  between  boundary  points 
and  points  not  belonging  to  the  mesh  are  to  be  set  equal 
to  zero.   We  now  transform  the  sum  through  partial  summa- 
tion, I.e.  we  arrange  the  sum  according  to  v,  and  split 
It  up  Into  a  siun  over  the  set  of  Interior  points,  G*  and 
a  sum  over  the  set  of  boundary  points,  ("l  .   Thus  we 
obtain: 

(1)  ^^11     B(u,v)  »  -h^  I  I  vL(u)  -  h  ^  vR(u). 

\  ^h  fh 

L(u)  Is  a  linear  difference  expression  of  second  order 
defined  for  all  Interior  points  of  C^: 

L(u)  »(au^)-  +  (bu^)7  +  (cUy)^  +  (dUy)y 

-  a\x^   -  pUy  +  (yu)-  +  (5u)y  -  gu  . 

R(u)  Is,  for  every  boundary  point,  a  linear  difference 
expression  whose  exact  form  will  not  be  given  here. 
If  we  arrange  5"  J!  B(u,v)  according  to  u,  we 
find  ^ 

(2)  h^  I  I  B(u,v)  =  -h^  ^  S   uM(v)  -  hX  uS(v). 

^h  «A  TA 

M(v)  is  called  the  adjoint  difference  expression  of  L(u) 
and  has  the  form 
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M(v)  =  (av^)-  +  (bVy)-  +  (cv^)-  +  (dvy)- 
+  (av)-  +  (Pv)-  -  Yv^  -  6Vy  -  gv. 


while  S(v)   is  a  difference  expression  corresponding  to 
R(u)   for  the  boundary. 

Formulas  (l)  and  (2)  give 

(5)   h  2!  Z  (vL(n)  -  uM(v))  +  h  ^   (vR(u)  -  uS(v))  =  0  . 

Formulas  (l),  (2),  and  (5)  are  called  Green's  Formula. 

The  simplest  and  most  Important  case  results  if 
the  bilinear  form  is  symmetric,  i.e.  if  the  relations 
b=c,  a«Y,  P=6  hold.   In  this  case  L(u)   is 
identical  with  M(u)  —  the  self-adJolnt  case  —  and  it 
can  be  derived  from  the  quadratic  expression 

2  2  2 

B(u,u)  =  au„  +  2bu  u  +  du   +  2au  u  +  2pu  u  +  gu  . 

In  the  following  we  shall  limit  ourselves  mainly  to 
expressions  L(u)  which  are  self-adjoint.   The  character 
of  the  difference  expression  L(u)   depends  principally 
on  the  nature  of  those  terms  in  the  quadratic  form  B(u,v) 
which  are  quadratic  in  the  first  difference  quotients.   We 
call  this  part  of  B(u,u)   the  characteristic  form; 

P(u,u)  «  au^^  +  2bu  u^  +  du^^  . 
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We  call  the  corresponding  difference  expression  L(u) 

elliptic  or  hyperbolic  depending  on  whether  the  function 

P(u,u)  of  the  difference  quotients  is  (positive)  definite 

or  Indefinite. 

The  difference  expression  Au  =  u  —  +  u  —  with  which 

XX    yy 

we  shall  concern  ourselves  in  the  following  paragraph  is 
elliptic.  I.e.  it  comes  from  the  quadratic  expression 

B(u,u)  =*  u^  +  Uy^  or  u-^  -(  Uy^  . 


The  corresponding  Green's  formulas  are 

(4)       h^  1    1   (u/  +  u  2)  =  -h^  ^  I  uAu  -  h  I  uR(u) 
%  ""  %  h 


(5)       h^  I  I  (vAu  -  uAv)  +  h  X  (vR(u)  -  uR(v))  =  0- 
^h  h 


3 


The  difference  expression,  Au  =  u  —  +  u  —  ,  Is  ob- 

XX    yy       p     _ 

vlously  the  analogue  of  the  differential  expression  *_H  +  oL^ 

^^    ^ 
for  a  function  u(x,y)  of  the  continuous  variables 

X  and  y.   Written  out  explicitly  the  difference 

expression  is 

Au  =  -^  Ju(x+h,y)+u(y,y+h)+u(x-h,y)+u(x,y-h)-4u(x,y)> 


^The  boundary  expression  R(u)  may  be  written  as  follows: 
Let  u  ,u,,...,Uy  be  values  of  the  function  at  a  boundary 

point  and  at  its  v  neighboring  points  (^0)>  then 

R(u)  =  ^  (uj+...+u^  -  v^Uq). 
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Therefore  "p"  ^^  is  the  excess  of  the  arithmetic  mean 
of  the  functional  values  at  the  four  neighborhood  points 
over  the  functional  value  at  the  point  in  question. 

Completely  similar  considerations  lead  to  linear 
difference  expressions  of  the  fourth  order  and  cor- 
responding Green's  formula,  provided  one  starts  from 
bilinear  difference  expressions  which  are  formed  from 
the  difference  quotients  of  second  order.   Consider  for 
example  the  difference  expression 


/\u  =  u  --  +  2u  -  -  +  u 

^-^    xxxx    xxyy   yyyy 


This  corresponds  to  the  quadratic  expression 


provided  one  orders  the  sum 


(J'h 


h 

according  to  v,   or  equlvalently  replaces  u  by  Au 
In  equation  (5).   One  must  notice  however  that  in  the 
expression  ^  /\u.   the  functional  value  at  a  point  Is 
connected  with  the  values  at  its  neighboring  points  and 
at  their  neighboring  points,  and  accordingly  is  defined 
only  for  such  points  of  the  region  G,   as  are  also  in- 
terior points  of  the  region  ^\*      (See  Section  5)»   The 
entirety  of  such  points  we  designate  as  G* '  ,   We  obtain 
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then  Green's  formxila 


(6)        h^Y"  Y_  A^-Av  =  h^  ^v.A  A^  +  li     S_  v.R(u)      , 

where  R  is  a  definable  linear  difference  expression 
for  each  point  of  the  boundary  strip   (7+  \,     ,       p, 
indicates  the  set  of  boundary  points  of  G'  . 

Section  2.  Boundary  Value  and  Eigenvalue  Problems 
1»   The  theory  of  boundary  value  problems 
The  boundary  value  problem  for  linear  elliptic 
homogeneous  difference  equations  of  second  order,  which 
corresponds  to  the  classical  boxindary  value  problem  for 
partial  differential  equations,  can  be  formulated  In 
the  following  way. 

Let  there  be  given  a  self-adjoint  elliptic  linear 
difference  expression  of  second  order,  L(u),  in  a  mesh 
region,  G.  ,  L(u)  results  from  a  quadratic  expression 
B(u,v)  which  is  positive  definite  in  the  sense  that  it 
cannot  vanish  If  u   and  u   do  not  themselves  vanish. 

A  fvmction,  u,   is  to  be  determined  satisfying 
in  G,   the  difference  equation 

L(u)  =  0   , 

and  assuming  prescribed  values  at  the  boxindary  points. 

Under  these  requirements  there  will  be  exactly 
as  many  linear  equations  as  there  are  interior  points 
of  the  mesh  at  which  the  function  u  Is  to  be 
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determined  ^.   Some  of  these  equations  which  Involve 
only  mesh  points  whose  neighbors  also  lie  In  the  In- 
terior of  the  region  are  homogeneous;   others  which  in- 
volve boundary  points  of  the  mesh  region  are  inhomo- 
geneous.   If  the  right-hand  side  of  the  inhomogeneous 
equations  is  set  equal  to  zero,  that  is  if  u  =  0  on 
the  bovindary,  then  it  follows  from  Green's  formula  (1), 
by  setting  u  =  v  that  B(u,u)  vanishes,  and  further 
from  the  definlteness  of  B(u,u)   that  u^  and  u^ 
vanish,  and  hence  that  u  itself  vanishes.   Thus  the 
difference  equation  for  zero  boundary  value?;  has  the 
solution  u  =  0,  or  in  other  words  the  solution  is 
uniquely  determined  since  the  difference  of  two  solu- 
tions with  the  same  bovindary  value  must  vanish.  Fur- 
ther, if  a  linear  system  of  equations  with  as  many  un- 
knowns as  equations  is  such  that  for  vanishing  right- 
hand  side  the  unknowns  must  vanish,  then  the  fundamental 
theorem  of  the  theory  of  equations  asserts  that  for  an 
arbitrary  right-hand  side  exactly  one  solution  must  ex- 
ist. In  our   case  there  follows  at  once  the  existence 
of  a  solution  of  the  botmdary  value  problem. 

Therefore  we  see  that  for  elliptic  difference 

^If  the  matrix  of  the  linear  system  of  equations  corres- 
ponding to  an  arbitrary  difference  equation  of  second  or- 
der,  L(u)=0,  is  transposed,  then  the  transposed  set  of 
equations  corresponds  to  the  adjoint  difference  equation 
M(v)=0.   Thus  the  above  self -ad joint  system  gives  rise  to 
a  set  of  linear  equations  with  symmetric  coefficients. 
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equations  the  uniqueness  and  existence  of  the  solution 
of  the  boundary  value  problem  are  related  to  each  other 
through  the  fundamental  theorem  of  the  theory  of  linear 
equations,  whereas  for  partial  differential  equations 
both  facts  must  be  proved  by  quite  different  methods. 
The  basis  for  this  difficulty  In  the  latter  case  is  to 
be  found  in  the  fact  that  the  differential  equations 
are  no  longer  equivalent  to  a  finite  number  of  equations, 
and  so  one  can  no  longer  depend  on  the  equality  of  the 
nxoraber  of  unknovms  and  equations. 

Since  the  difference  expression 

can  be  derived  from  the  positive  definite  quadratic 
expression 

2  r-  T-  /..2  ^  „2. 


h'^^(u^.u^ 


) 


the  boiindary  value  problem  for  the  difference  expres- 
sion is  uniquely  solvable. 

The  theory  for  difference  equations  of  higher 
order  is  developed  in  exactly  the  same  way  as  that  for 
difference  equations  of  second  order;  for  example  one 
can  treat  the  fourth  order  difference  equation 

A  A  ^  =  0    • 
In  this  case  the  values  of  the  fimction  must  be  prescribed 
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on  the  boundary  strip  [^,    +  \,     ,      Evidently  here  also 
the  difference  equation  yields  just  as  many  linear  equa- 
tions as  there  are  unknown  functional  values  at  the 
points  of  G^'il*   In  order  to  demonstrate  the  existence 
and  uniqueness  of  the  solution  one  needs  only  to  show 
that  a  solution  which  has  the  value  zero  in  the  bound- 
ary strip  \~^,    +    Pj^  necessarily  vanishes  identically. 
To  this  end  we  note  that  the  sum  over  the  corresponding 
quadratic  expression 

(7)  ^^ZZ  ^^""^^ 

^"h 

for  such  a  function  vanishes,  as  can  be  seen  by  trans- 
forming the  sum  according  to  Green's  formula  (6).  The 
vanishing  of  the  sura  (?)   however  implies  that  /\u 
vanishes  at  all  points  of  G'  ,   and  according  to  the 
above  proof  this  can  only  happen  for  vanishing  bound- 
ary values  if  the  fvmction  u  assumes  the  value  zero 
throughout  the  region.   Thus  our  assertion  is  proved, 
and  both  the  existence  and  uniqueness  of  the  solution 
to  the  boundary  value  problem  for  the  difference  equa- 
tions are  guaranteed^. 


■P- — . 

"^For  the  actual  process  of  carrying  through  the  solu- 
tion of  the  boundary  value  problem  by  an  iterative 
method,  see  among  others  the  treatment:  "Uber  Randwert- 
anjgahen  bei  partieller  Dlfferenzenglerthunger"  by  R. 
Courant,  Zeltschr.  f.  angew,   Mathematik  u,  Mechanik  6 
(1925),  pp.  32^-3^5'      Also  there  is  a  report  by  H.  Henky, 
in  Zeitschr,  f.  angew.  Math.  u.  Mech.  2  (1922),  p.  58  ff. 
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2.   Relation  to  the  mlnimmn  problem 
The  above  boiindary  value  problem  is  related  to 
the  following  minimum  problem:  among  all  functions 
i|?(x,y)   defined  in  the  mesh  region  G,   and  assuming 
given  values  at  the  boundary  points,  that  function 
^  =  u(x,y)   is  to  be  found  for  which  the  sum 


BiLff  ,if) 


\ 


over  the  mesh  region  assumes  tho  least  possible  value. 
We  assume  that  the  qviadratic  difference  expression  of 
first  order,  B(u,u)  is  positive  definite  in  the  sense 
described  in  Section  1,  part  2.   One  can  show  that  the 
difference  equation  L(  (p  )  =  0  results  from  this  mini- 
mum requirement  on  the  solution  (f  -   u(x,y),  where 
L{^  )      Is  the  difference  expression  of  second  order  de- 
rived previously  from  B((^,i^).   In  fact  this  can  be 
seen  either  by  applying  the  r\iles  of  differential  cal- 
culus to  the  s\ims  h*^  V"  B(Cf?,<:^)  as  a  function  of 

\ 

a  finite  number  of  values  of  cp    at  the  grid  points,  or 
by  employing  the  usual  methods  from  the  calculus  of 
variations. 

By  way  of  example,  solving  the  boundary  value 
problem  of  finding  the  solution  to  the  equation  ^(^  =  0 

which  assumes  given  boxuidary  values,  is  equivalent  to 

o   2     2 

minimizing  the  sum  h  ^  y"(^  ^  •*■  ^     )      over  all  functions 
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assTOBilng  given  values  on  the  boundary  strip   iL*   Be- 
sides this  sum  there  ara  yet  other  quadratic  expres- 
sions in  the  second  derivatives  which  give  rise  to  the 
equation  ^  /\u.  =   0  under  the  process  of  being  mini- 
mized. For  example  this  Is  true  in  G,   for  the  sum 

That  the  minimum  problem  posed  above  always  has 
a  solution  follows  from  the  theorem  that  a  continuous 
function  of  a  finite  n\imber  of  variables  (the  function- 
al values  of  ^f  at  the  grid  points)  always  has  a  mini- 
mum if  it  is  bounded  from  below  and  if  it  tends  to  in- 
finity as  soon  as  any  of  the  Independent  variables  goes 
to  infinity^. 

3»   Green's  function 

It  is  possible  to  treat  the  boundary  value  prob- 
lem for  the  Inhomogeneous  equation,  L(u)  =  -f,   in 
much  the  same  way  as  the  homogeneous  case,   L(u)  =  0, 
For  the  Inhomogeneous  case  it  is  sufficient  to  consider 
only  the  case  of  zero  boundary  conditions,  since  dif- 
ferent boundary  conditions  can  be  taken  care  of  by  ad- 
ding a  suitable  solution  of  the  homogeneous  equation. 
To  solve  the  linear  system  of  equations  representing 

It  can  easily  be  verified  that  the  hypotheses  for  the 
application  of  this  theorem  are  satisfied. 

-  15  - 


NYO-7689 

the  boundary  value  problem,   L(u)  =  -f,   we  first  choose 
as  the  function  f(x,y)   a  function  which  assuines  the 

value ^  at  the  point  x  =  4,  y  =  y  of  the  mesh. 

If  K(x,y{4,'n)   is  the  solution  (vanishing  on  the 
boundary)  of  this  difference  equation  which  depends  on 
the  parametric  point  (^,^)t      then  the  solution  for  an 
arbitrary  boundary  condition  can  be  represented  by  the 
sum 

u(x,y)  =  ^^     y  ^^^>'^»^»\  ^   ^^^»^  )   • 

Tlie  function  K(x,y;4,'>^  )  which  depends  on  the 
points   (x,y)   and   (C,  >i )   is  called  the  Green's  func- 
tion of  the  differential  expressions  L(u).   If  we  call 
the  Green's  function  for  the  adjoint  expression  M(v), 
K(x,y ;£;,"*]  ),   then  the  equivalence 

Kd,")!  54, n  )  =  K(C,">i  il,\  )         holds, 
as  can  be  seen  to  follow  immediately  from  Green's  form- 
ula (3)  when  u  =  K(x,yj4,'>2  ^»  ^^^     ^  ~   ^^^»7»^*^  )• 
For  a  self-adJolnt  difference  expression  the  above  re- 
lation gives  the  symmetric  expression: 

if.  Eigenvalue  problems 

Self-adjoint  difference  expressions,  L(u),  give 
rise  to  eigenvalue  problems  of  the  following  type:   find 
the  value  of  a  parameter  X,   the  eigenvalue,  such  that 
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in  G,  ,   a  solution,  the  eigenfu  notion  can  be  found 
for  the  difference  equation 

L(u)  +  Xu  =  0   , 

where  u  Is  to  be  zero  on  the  boundary,  \-^' 

The  eigenvalue  problem  is  equivalent  to  finding 
the  principle  axes  of  the  quadratic  form  B(u,u). 
Exactly  as  many  eigenvalues  and  corresponding  eigen- 
functions  exist  as  there  are  interior  mesh  points  of 
the  region  G,  .   The  system  of  eigenfunctions  and  eigen- 
values (and  a  proof  of  their  existence)  is  given  by  the 
minimum  problem: 

Among  all  functions,  (f>(x,y),  vanishing  on  the 
boundary,  and  satisfying  the  orthogonality  relation 

L^^^  =  0  ,   (v=l,...,m-l) 
'h 
and  normalized  such  that 


hV  ^  cp  u' 


XL^'-^ 


h' 
the   function,     Cp  =  u,      is  to  be   found  for  which  the    svun 


h^  y  Bicp  ,cp) 


^h 

assumes    its  minimum  value.      The   value   of  this  miniraxun 
is   the     m         eigenvalue,    and  the  function  for  which  it 
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1"  Vi  7 

is  e.ssuined   is   the     m         eigenfunction   , 


7 
From  the   orthogonality  condition  on  the  eigenfimctlons. 


h^y  ^  u^^^u^^^)  =  0    ,    (v^; 


it  follows  that  each  function,  g(x,y),  which  vanishes 
on  the  boundary  of  the  mesh  can  be  expanded  in  terras  of 
the  eigenfunctions  in  the  form 

v=l 

where  the  coefficients  are  determined  from  the  equations 


\^ 


In  this  way  in  particular  the  following  represent- 
ation for  the  Green's  functions  maY  be  derived. 


^   ,(v).„..(v) 
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o 

Section  3»   Connections  with  the  Problem  of  the  Random  Walk 

The  theme  of  the  following  Is  related  to  a  ques- 
tion from  the  theory  of  probability,  namely  the  problem 

9 
of  the  random  walk  In  a  bounded  region.   We  consider 

the  lines  of  a  mesh  region  G.   as  paths  along  which  a 

particle  can  move  from  one  grid  point  to  a  neighboring 

one.   In  this  net  of  streets  the  particle  can  wander 

aimlessly,  and  it  can  choose  at  random  one  of  the  fo\ir 

directions  leading  from  each  intersection  of  paths  of 

the  net  -  all  four  directions  being  equally  probable. 

The  walk  ends  as  soon  as  a  boxindary  point  of  G,   is 

reached  because  here  the  particle  must  be  absorbed. 

We  ask: 

1)  What  is  the  probability  W(PjR)  that  a  ran- 
dom walk  starting  from  a  point  P  reaches  a 
particle  point  R  of  the  boundary? 

2)  What  is  the  mathematical  expectation  v(P;Q,) 
that  a  random  walk  starting  from  P  reaches 
a  point  Q,  of  G,   without  touching  the 
boundary? 

Q 

Section  3  is  not  prerequisite  to  Section  ij.. 

9 
The  present  treatment  Is  essentially  different  from  the 

familiar  treatments  which  can  be  carried  through,  say 

for  example  in  the  case  of  Brownian  motion  for  molecules. 

The  difference  lies  precisely  In  the  way  in  which  the 

boundary  of  the  region  enters. 
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This  probability  or  mathematical  expectation,  re- 
spectively will  be  defined  more  precisely  by  the  follow- 
ing process.  Assume  that  at  the  point  P  there  is  a 
unit  concentration  of  matter.   Let  this  matter  diffuse 
into  the  mesh  with  constant  velocity,  traveling  say  a 
mesh  width  in  unit  time.  At  each  meshpoint  let  exactly 
one  fourth  of  the  matter  at  the  point  diffuse  outwards 
in  each  of  the  four  possible  directions.   The  matter 
which  reaches  a  boundary  point  is  to  remain  at  that 
point.   If  the  point  of  origin  P  is  itself  a  boundary 
point,  then  the  matter  never  leaves  that  point. 

We  define  the  probability  w(P,R)   that  a  ran- 
dom walk  starting  from  P  reaches  the  boundary  point 
R  (without  previously  attaining  the  boundary),  as  the 
amount  of  matter  which  accumulates  at  this  boundary 
point  over  an  infinite  amount  of  time. 

We  define  the  probability  E  (PjQ,)   that  the  walk 
starting  from  the  point  P  reaches  the  point  Q  In 
exactly  n  steps  without  touching  the  boundary  by  the 
amount  of  matter  which  accumulate  In  n  vmits  of  time 
provided  P  and  Q,  are  both  Interior  points.   If 
either  P  or  Q  are  boundary  points  then  E   is  set 
equal  to  zero. 

The  value  E^(P;Q)   is  exactly  equal  to  lA"^  times 
the  number  of  paths  of  n  steps  leading  from  P  to  Q 
without  touching  the  bovindary.   Thus  E  (PjQ)  =  E  (Q;P). 
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We  define  the  mathematical  expectation  V(P;Q) 
that  one  of  the  paths  considered  above  leads  from  P 
to  the  point  Q,  by  the  infinite  svun  of  all  of  these 

possibilities, 

00  10 

V(P;Q)  =y  E^(P;Q) 

i.e.,  for  the  interior  points  P  and  Q,   the  svun  of 
all  the  concentrations  which  have  passed  through  the 
point  Q  at  different  times.   This  will  be  assigned 
the  expected  value   1  for  a  concentration  originating 
at  Q. 

If  the  amount  arriving  at  the  boundary  point  R 
in  exactly  n  steps  is  designated  as  P  (PjR),  then 
the  probability  w'P;R)   is  given  by  the  series 

00 

w(PjR)  =  y  P^(FiR)   . 
v=0 

All  the  terms  of  this  series  are  positive  and 
the  partial  svim  is  bounded  by  one  (since  the  concentra- 
tion reaching  the  boundary  can  be  made  up  of  only  part 
of  the  initial  concentration),  and  therefore  the  con- 
vergence of  the  series  is  assured* 

Now  it  Is  easy  to  see  that  the  probability  E  (P;Q), 
that  is,  the  concentration  reaching  the  point  Q  In 
exactly  n  steps  tends  to  zero  as  n  increases.  For 
if  at  any  point  Q,  from  which  the  boundary  point  R 

The  convergence  will  be  shown  below. 
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can  be  reached  In  m  steps,  we  have  E  (PjQ)  >  a  >  0, 
then  at  least  a/\^  of  the  concentration  at  Q  will 
reach  the  point  R  after  in  steps.   However,  since 

00 

the  sum  T^  F  (PjR)   converges,  the  concentration  reach- 

v=0  ^ 
ing  R   goes  to  zero  with  increasing  time,  where  the 

value  of  E  (PjQ)  must  itself  vanish  as  time  increases: 
n 

that  is  the  probability  of  an  infinitely  long  walk  re- 
maining in  the  Interior  of  the  region  is  zero. 

Prom  this  it  follows  that  the  entire  concentra- 
tion eventually  reaches  the  boundary;  or  in  other  words 

that  the  siim  ^  w(PjR)   over  all  the  boundary  points 

R 
R  is  equal  to  one. 

The  convergence   of  the   infinite   series  for  the 

mathematical  expectation     v(P;Q) 

OD 

v(P;Q)  =  y  E^(P;Q) 

remains  to  be  shown. 

To  this  end  we  remark  that  the  quantity  E  (PjQ) 
satisfies  the  following  relations 

where  Q,   through  Q,   are  the  four  neighboring  points 

of  the  Interior  point  Q,.   That  is,  the  concentration 

st  1 

at  the  point  Q  at  the   n+1    step  consists  of  y 

times  the  s um  of  the  concentrations  at  the  four  neigh- 

boring  points  at  the   n    step.   If  one  of  the  neighbors 
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of  Q,   for  example   ^-=  R,   is  a  boundary  point  then  it 

follows  that  no  concentration  flows  from  this  boundary 

point  to  Q   since  the  expression  E  (P;R)   is  zero  in 

this  case.   Furthermore,  for  an  interior  point,   E^(PjP)=l 

and  of  course   E^(PjQ)=0. 

Prom  these  relationships  we  find  for  the  partial  sum 

v^(PiQ)  =  Y^   2^(PiQ) 
v=0 
the  equation 

^n^l^P'^)  =t[v^'^l)-^V^'^2)^V^iS^-*-V^'\)  ]   ' 
for  P  ^  Q  .   For  P  =  Q, 

v^^^CPjP)  =  1^  w^(PjP^)-Kv^(P;P2)H-v^(P;P3)+v^(P;P^)  j   , 

that  is,  the  expectation  that  a  point  goes  back  into  it- 
self consists  of  the  expectation  that  a  non-vanishing 
path  leads  from  P  back  again  to  itself  -  namely 
^  |v^(P;P^)+v^(PjP2)+v^(P;P3)+v^(P;P^)  |  ,  together 

with  the  expectation  unity  that  expresses  the  initial 
position  of  the  concentration  itself  at  this  point. 

The  quantity  v  (PjQ)   thus  satisfies  the  follow- 
ing difference  equation   . 

A^n^Pj^)  =  ^   E^(P;^)  »   f or  P  ?^  Q  , 
h 

Av„(P;Q)  =  -^CE^CPjQ)-!)  ,   for  P  =  Q  . 
n       j^^  n 

This  defines  the  /Vope ration  for  the  variable  point  Q. 

This  equation  can  be  interpreted  as  an  equation  of  the 
heat  conduction  type.   That  is  if  the  function  v  (P;Q) 
is  considered,  not  as  a  function  of  the  index  n  as  in 
our  presentation  above,  but  rather  as  a  function  of  time, 
t,  which  is  proportional  to  n,   so  that  t  =  nT   and 
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V  (P;Q)  is  equal  to  zero  when  Q   is  a  boundary  point. 

The  solution  of  this  bovmdary  value  problem  for  arbi- 
trary right-hand  side  is  viniquely  determined  as  we  have  ex- 
plained earlier  (Section  2,    Part  1.),  and  depends  continu- 
ously on  the  right-hand  side.   Since  the  variables  E  (PjQ) 

tend  to  zero,  the  solution  v  (P:Q)   converges  to  the  solu- 
'  n 

tion  v(PjQ,)   of  the  difference  equation 
Av(P;Q)  =  0    for   P  ^  Q 
Av(PjQ)  =  -  -^  for  P  =  Q, 

with  boundary  values  v(P;R)  =  0, 

Thus  we  see  that  the  mathematical  expectation  v(P;Q) 
exists  and  is  non  other  than  the  Green's  function  for  the 
difference  equation  /\u  =  0»   except  for  a  factor  of  \\., 
The  symmetry  of  the  Green's  fxinction  is  an  immediate  conse- 
quence of  the  symmetry  of  the  quantity  E  (PjQ)  which  was 
used  to  define  it. 

The  probability  w(P;R)   satisfies,  with  respect  to 
P,  the  relation 

w(P;R)  =  -^  U(P3_;R)+w(P2;R)+w(P2;R)+w(P^;R)  j  , 
and  thus  the  difference  equation 

Aw  =  0 

V  (P;Q)  =  v(P;Qjt)  =  v(t),   then  the  above  equations  can  be 
written  in  the  following  fom: 

^v(t)  =  iq  .  '''^^V'"'  f-  p  ^ «  . 

h^        ^ 

^wt)  .^  .^v(t.r_)-v(t)  .i)     ,„^   ,^^    _ 

For  a  similar  difference  equation  which  has  a  para- 
bolic differential  equation  as  its  limiting  form,  see  Section 
6  of  the  second  half  of  the  paper. 
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That  Is,  if   P^,P2,P^,P.  ,   are  the  four  neighboring 
points  of  the  interior  point   P,   then  each  path  from 
P   to  R  must  pass  through  one  of  these  four  direc- 
tions, and  each  of  the  four  is  equally  likely.   Further- 
more, the  probability  of  going  from  one  boundary  point 
R  to  another  R'   is  w(RjR')  =  0  xinless  the  two 
points  R  and  R'   coincide,  in  which  case  w(RjR)  =  1. 
Thus  w(PjR)   is  that  solution  of  the  boundary  value 
problem  /\w  =  0  which  assumes  the  value  1  at  the 
boundary  point  R  and  the  value   0  at  all  other 
points  of  the  boundary.   Therefore  the  solution  of  the 
boundary  value  problem  for  an  arbitrary  boundary  value 

u(R)   has  the  simple  form  u(P)  =  Yl   w(P;R)u(R),  wher<^ 

R  ^2 

the  sum  is  to  be  extended  over  all  the  boundary  points   . 

If  the  function  u  =  1  is  substituted  for  u  in  this 

expression,  then  we  check  the  relation  1  =  1^  w(PjR)« 

R 
The  interpretation  given  above  for  Green's  func- 
tion as  an  expectation  yields  immediately  further  prop- 
erties. We  mention  only  the  fact  that  the  Green's  func- 
tion decreases  if  one  goes  from  the  region  G  to  a 
subregion  "3  lying  within  G;   that  is  the  number  of 
possible  paths  for  steps  on  the  mesh  leading  from  one 

12  ' 

Moreover  it  is  easy  to  show  that  the  probability  w(PjR) 

of  reaching  the  boundary  is  the  boundary  expression 

3i»(K(P,Q),  constructed  from  the  Green's  function  K(P;Q)  In 

terms  of  Q,  where  u(x,y)  is  to  be  identified  with  w(P,Q), 

and  v(x,y)  with  v(P,Q)  in  Green's  formula  (5). 
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point  P  to  another  Q  (without  touching  the  bovindary), 
decreases  as  the  region  decreases. 

Of  course  corresponding  relationships  hold  for 
more  than  two  Independent  variables.  We  note  only  that 
other  elliptic  difference  equations  admit  a  similar 
probability  interpretation. 

If  the  limit  for  vanishing  mesh  width  is  considered 
by  methods  given  in  the  following  section,  then  the  Green's 
fiinction  on  the  mesh  goes  over  to  the  Green's  function 
of  the  potential  equation  except  for  a  numerical  factor; 

a  similar  relationship  holds  between  the  expression 

w( P*R ) 

,'  •'-     and  the  normal  derivative  of  the  Green's  func- 
tion at  the  boundary  of  the  region.   In  this  way  for  ex- 
ample the  Green's  function  for  the  potential  equation 
could  be  Interpreted  as  the  specific  mathematical  ex- 
pectatlon  of  going  from  one  point  to  another,  "^  without 
reaching  the  bo^indary. 

In  going  over  to  the  limit  of  a  continuum  from 
the  mesh,  the  influence  of  the  direction  in  the  mesh  pre- 
scribed for  the  random  walk  vanishes.   This  fact  suggests 
that  it  would  be  of  some  Interest  to  consider  limiting 
cases  of  more  general  random  walks  for  which  the  limit- 
ations on  the  direction  of  the  walk  are  relaxed.  This 
lies  outside  of  the  scope  of  this  presentation  however 
and  we  can  only  hope  to  renew  the  question  at  some  other 
opportiinity. 


''Here  the  a  priori  expectation  of  reaching  a  certain  area 
element  is'' understood  to  be  equal  to  the  area  of  the  element, 
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Section  k»      The  Solution  of  the  Differential  Equation  as 
a  Limiting  Form  of  the  Solution  of  the  Difference  Equation 

1.   The  boundary  value  problem  of  potential  theory. 

In  letting  the  solution  of  the  difference  equation 
tend  to  the  solution  of  the  corresponding  differential 
equation  we  shall  relinquish  the  greatest  possible  degree 
of  generality  with  regard  to  the  boundary  and  boxindary 
values  in  order  to  demonstrate  more  clearly  the  charac- 
ter of  our  method  ^.   Therefore  we  assume  that  we  are 
given  a  simply  connected  region  G     with  a  boundary  formed 
of  a  finite  number  of  arcs  with  continuously  turning  tan- 
gents.  Let   f(x,y)   be  a  given  function  which  is  contin- 
uous and  has  continuous  partial  derivatives  of  first  and 
second  order  in  a  region  containing   G.   If  G,   Is  the 
mesh  region  with  mesh  width  h  belonging  to  the  region 
G,      then  let  the  boundary  value  problem  for  the  difference 
equation  /\u  =  0  be  solved  for  the  same  boundary  values 
which  the  function  f(x,y)   assumes  on  the  boundary;   let 
u,  (x,y)   be  the  solutlonc.   We  shall  prove  that  as  h — >0 
the  function  u,  (x,y)   converges  to  a  function  u(x,y) 
which  satisfies  in   G  the  partial  differential  equation 

^^^—p   "*"  ''^^  ~  ^  ^^'^  takes  the  value  of  f(x,y)   at  each 

<^x    <)y 

of  the  points  of  the  boundary.   We  shall  show  further  that 

^We  mention  however  that  carrying  through  our  method  for 
more  general  boundaries  and  boundary  values  in  no  way 
causes  any  fundamental  difficulty. 
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for  any  region  lying  entirely  within   G  the  difference 
quotients  of  u,   of  arbitrary  order  tend  uniformly  to- 
wards the  corresponding  partial  derivatives  of  u(x,y). 
In  the  convergence  proof  it  is  convenient  to  re- 
place the  boundary  condition  u  =  f  by  the  weaker  re- 
quirement that 

(u-f)^dxdy  -->  0  as  r  — >  0   , 


where   S  is  that  strip  of  G  whose  points  are  at  a 
distance  from  the  boundary  smaller  than  r. 


The  weaker  boundary  value  requirement  does  in  fact 
provide  the  unique  characterization  of  the  solution, 
as  can  be  seen  from  the  easily  proved  theorem:   If 
the  boundary  condition  above  is  satisfied  for  f(x,y)=G 
for  a  function  satisfying  the  equation 


in  the  interior  of   G  and  if      Ht^I      +   (4^)   /dxdy 


exists,  then  u(x,y)  is  identically  zero.  (See  Courant, 
"Uber  die  Losung  der  Differentialgleichungen  der  Physik", 
Math  Ann.  8^  exp.  p.  296  ff . ) 

In  the  case  of  two  independent  variables  the  bovuid- 
ary  values  are  actually  attained,  as  follows  from  the 
weaker  requirement;   but  in  the  case  of  more  variables 
the  corresponding  result  cannot  in  general  be  expected 
since  there  can  exist  exceptional  points  on  the  bound- 
ary at  which  the  boundary  value  is  not  taken  on  even 
though  a  solution  exists  under  the  weaker  requirement. 
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The  convergence  proof  depends  on  the  fact  that  for  any 
subreglon  G-   lying  entirely  within  G,   the  function 
u,  (x,y)   and  each  of  its  difference  quotients  is  bound- 
ed and  tiniformly  continuous  as  h  -->  0   in  the  follow- 
ing sense:  For  each  of  these  functions,  say  w,  (x,y), 
there  exists  a  5(e)   depending  only  on  the  subregion 
and  not  on  h  such  that 

|wj^(P)  -  w^(P3^)|  <  e 

provided  the  mesh  points  P  and  P,  lie  In  the  same 
subregion  of  G,  and  are  separated  from  each  other  by 
a  distance  leas  than  5(e). 

Once  uniform  continuity  In  this  sense  (equl- 
contlnulty)  has  been  established  we  can  in  the  usual 
way  select  from  the  functions  u,   a  subsequence  of 
functions  which  tend  uniformly  In  any  subregion  G 
to  a  limit  function  u(x,y),  while  the  difference  quo- 
tients of  u,   tend  uniformly  towards  the  partial  deriv- 
atives of  u»   The  limit  ftinctlon  then  possesses  deriv- 
atives of  arbitrarily  high  order  in  any  proper  subregion 

*  2 

G   of  G  and  satisfies  ^  u  =  0  In  this  region.   If 

we  can  show  also  that  u  satisfies  the  boundary  condi- 
tion we  can  regard  It  as  the  solution  of  our  boundary 
value  problem  for  the  region  G.   Since  this  solution 
la  uniquely  determined.  It  is  clear  that  not  only  a  par- 
tial sequence  of  the  functions  u,  ,  but  this  sequence 
of  functions  Itself  possesses  the  required  convergence 
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properties. 

The  uniform  continuity  (equi-contlnuity )  of  our 
quantities  may  be  established  by  proving  the  following 

lemmas. 

1)  As  h  — >  0  the  sums  over  the  mesh  region 

h^  E  E  "^  and  h^  E  E(^x  *   "t' 

remain  bounded. 

2)  If     w  =  w,       satisfies   the   difference   equa- 
tion    Aw  =0     at  a  mesh  point    of     G^,      and 
if,   as     h  — >  0     tiie    svm. 

2  T-  r—  .  2 


'^  LI" 


.-> 


extended  over  a  mesh  region  G^  associated 
with  a  subregion  G   of  G  remains  bounded, 
then  for  any  fixed  subregion  G    lying  en- 
tirely within  G   the  sum 


^'  I  I'"x  ^  "y 


4Ht 


.-rc-$* 


over  the  mesh  region  G'  associated  with 
G**  likewise  remains  bo\inded  as  h  — >  0. 
Using  (1)  there  follows  from  this,  since 
all  of  the  difference  quotients  w  of  the 


^^Here  and  in  the  following  we  drop  the  index  h  from 
the  grid  functionSf 
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function  u,   again  satisfy  the  difference 


equation  /\w  =  0,   that  each  of  the 


sums 


h  21  ]^  w   is  bounded. 

3)   From  the  boundedness  of  these  sums  there 

follows  finally  the  boundedness  and  uniform 
continuity  of  all  the  dlfft.rence  quotients 
themselves. 

2.   Proof  of  the  lemmas 

The  proof  of  1)   follows  from  the  fact  that  the 
fvmctlonal  values  u.   are  themselves  bounded.   For  the 

greatest  (or  least)  value  of  the  function  is  assumed  on 

17 
tiie  boundary    and  so  is  boxinded  by  a  prescribed  fin- 
ite value.   The  boundedness  of  the  sum  hy"  TZint  "*"  ^-) 
Is  an  Immediate  consequence  of  the  minimxim  property  of 
our  mesh  function  formulated  in  Part  2  of  Section  2 
which  gives  in  particular 

but  as  h  — >  0  the  sum  on  the  right  tends  to  the  integral 
G  ^ 

l7 — — 

We  note  however  with  a  view  to  carrying  over  the  method 

to  other  differential  equations,  that  we  can  relax  these 
conditions.   To  this  end  we  need  only  to  bring  into  play 
the  Inequality  (l5)  or  to  use  the  reasoning  of  the  altern- 
ative (see  Part  l^.,    Section  l^) , 
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which,   by  hypothesis,   exists. 

To  prove      2)      we   consider   the   quadratic   sum 


int.Q^ 

where  the  siunmation  extends  over  all  the  interior  points 
of  a  square  Q,,   (see  Pig.  1).  We  call  the  values  of 

the  fvmction  on  the  bound- 
ary S,   of  the  square  Q-, 


^n 

h 

^0 

w 


1' 

ary     S 


and  those   on  the  bound- 


^,      of     Q^,   w^.      Then 


Fig.   1 


Green's  formula  gives 


( 8 )        intV^^ 


jy^  -  ^w^-  Vw' 


where  S^  and  S^  are  respectively  the  boundary  of  Q- 
and  the  square  boundary  of  the  lattice  points  Q   lying 
within  S. ,  while  C.   consists  of  the  four  corner  points 
of  the  bovmdary  of  Q,  . 

We  now  consider  a  sequence  of  concentric  squares 
Qq,Q^,«..,Qj^  with  boundaries  S^,S^, . . ,  ,S^,   where  each. 
boundary  is  separated  from  the  next  by  a  mesh  width. 
Applying  the  formula  to  each  of  these  squares  and  ob- 
serving that  we  have  always 
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2h^  ^(w^  +  w|)  <  h^  ^(w|  +  w|  +  w|)  ,   (k  >  1) 


we  obtain 
2h' 


^  ^(w^  +  w|)  <  Vw^  -  XL  '  Z"^   (1  <  k  <  n] 

Q.  S,_     S,.  ,    C" 


^-  -1,     "ic-1   ''k 


We  strengthen  the  inequality  by  neglecting  the  last 

non  positive  term  on  the  right  and  we  then  add  the  n 
inequalities  to  obtain 


2nh' 


(1  s      s      s 

^o  n      o      n 


Svumuing  this   Inequality  from     n  =  1     to     n  =  N 
we   get 

T2^2r-  T—   /    2   ^     2- 


Diminishing  the  mesh  width  h  we  can  make  the 
squares  Q   and  Q   converge  towards  two  fixed  squares 
lying  within  G  suid  having  corresponding  sides  separ- 
ated by  a  distance  a.   In  this  process  Nh  — >  a  and 
we  have,  independent  of  the  mesh  width 

(9)        '^'i:  E'-x  *  "'>  ^  7  r  z "'  • 


N 
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For  sufficiently  small  h  this  Inequality  holds 
of  course  not  only  for  two  squares  Q   and  Q^  but 
with  a  change  in  the  constant,   a,   for  any  two  sub- 
regions  of  G  such  that  one  is  contained  entirely 
within  the  other.   Thus  lemma   (2)   is  proved. 

In  order  to  prove  the  third  result,  that  vu 
and  all  its  partial  difference  quotients  w,   remain 
botmded  and  xiniformly  continuous  as  h  — >  0,  we  con- 
sider a  rectangle  R  with  corners  P  ,Q  ,P,Q  and 

o  o 

with  sides   Pq'^q  and  PQ  of  length  a  parallel  to 
the  x-a_xis. 

We  start  with  the  relation 


w(Q^)  -  w(p^)  =  hyw^-h^j;  J^w^ 


1  fi 

If  we  do  not  assume  that  Aw  =0,   then  in  place  of 
the  inequality   (9)  we  find; 


(10)   h^  ^  (w^  +  w^)  <  c^h^  ^w^  +  C2hV  Y}A^)^ 
G*^  G*  G* 


for  suitable  constants  c,   and  c,  independent  of  h, 
where   G**  lies  entirely  within  G*  and  G*  in  turn 
Is  contained  in  the  interior  of  G. 
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(11)  |w(Qq)    -  w(Pq)|    <  h^    |w^ 

PQ 


^"l.  Z  I V 


o   o 


which  is  a  consequence  of  it. 

Ae    then  let  the  side   PQ   of  the  rectangle  vary 

between  an  initial  line   Pi^t »   ^  distance   b   from 

P^Q,^  and  a  final  line  ^p^Z     ^   distance  2b  from 

P  Q,  ,   and  we  sum  the  corres- 
o  o 

ponding  t:  +  1  inequalities 
(11).   We  obtain  the  estimate 


|w(P,)  -w(Q^)|  <-igh2^^  Iw^l 


Ro 


Fig.  2. 


where  the  summations  extend  over  the  entire  rectangle, 

R^  =  P  Q,  PoQ.o»  Prom  Sohwarz's  Inequality  It  then  fol- 
2     O  O  <i  <; 

lows  that, 


:i2)    Iw(P^)  -  w(Q^)|  <  ^/2Sb  /h^j-  Y^l   +  /25b  /h Y^  Yy 

Ro  R« 


xy 


Since,  by  hypothesis,  the  sums  which  occur  here 

2 
multiplied  by  h   remain  botinded,  it  follows  that  as 

a  -••>  0  the  difference   |w(Pq)  -  w(Q^)  |  -->  0  inde- 
pendently of  the  mesh-width,  since  for  each  subreglon 
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G   of  G  the  quantity  b  can  be  held  fixed.   Con- 
sequently the  vmlform  continuity  (equi-continuity)  of 
w  =  w.   Is  proved  for  the  x-direction.   Similarly  it 
holds  for  the  y-direction  and  so  also  for  any  sub- 
region  G   of  G,   The  boundedness  of  the  function 
w,   in  G   finally  follows  from  its  uniform  continu- 


ity (equi-continuity)  and  the  boundedness  of  h^"  ^wf 

By  this  proof  we  establish  the  existence  of  a 
subsequence  of  functions  u.   which  converge  towards 
a  limit  function  u(x,y)  and  which  do  so  uniformly 
together  with  all  their  difference  quotients,  in  the 
sense  discussed  above  for  every  interior  subregion  of 
G.   This  limit  function  u(x,y)   has  throughout  G  con- 
tinuous partial  derivatives  of  arbitrary  order,  and 
satisfies  there  the  potential  equation: 

4-4  =  0  . 

3.   The  boundary  condition 

In  order  to  prove  that  the  solution  satisfies  the 
boundary  condition  formulated  above  we  shall  first  of 
all  establish  the  inequality 

(13)     h^^  X]  ^^  ^  ^^^^^  Y}^^  "■  ""y^  *  ^""^  Z  ""^ 

^r,h  ^r,h  Th 

where   S  ^  is  that  part  of  the  mesh  region  G,   which 
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lies  inside  a  boundary  strip  S  ♦   This  boiindary  strip 
S   consists  of  all  points  of  G  whose  distance  from 
the  boundary  is  less  than  r;   it  is  bounded  by   p 

and  another  curve   P  .   The  constants  A  and  B  de- 

r 

pend  only  on  the  region  and  not  on  the  function  v 
nor  the  mesh  width  h. 

In  order  to  prove  the  above  Inequality,  we  divide 
the  boundary,   P,   of  G  into  a  finite  number  of 
pieces  for  which  the  angle  of  the  tangent  with  one  of 
the  X  or  y-axes  is  greater  than  some  positive  val- 
ue (say  30  )•   Let  y,   for  Instance,  be  a  piece  of 

r 

r 


Fig.  3. 


Pig.  i^.. 

P     which  Is   sufficiently  steep  (In  the  above   sense) 
relative   to  the  x-axis,    (See  Fig.   14.),     Lines  parallel 
to  the  x-axls   intersecting     y     will  cut  a  section     y 
from  the  neighboring  curve      P   ,     and  will  define   to- 
gether with     y    and     y^,     a  piece     s         of  the  boiindary 
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strip  S  .   We  use  the  symbol  s  ,   to  denote  the 
r  "  r>h 

portion  of  G,   contained  in  s   and  denote  the  as- 
h  r 

sociated  portion  of  the  bovmdary   ["*,   by  y,  . 

We  now  imagine  a  parallel  to  the  x-axis  to  be 

drawn  through  a  mesh  point  P,   of  s  ,  ,   Let  it 

meet  the  boundary  Yv.  in  a  point  P.  .   The  portion 

of  this  line  which  lies  in  s  ,   we  call  p  ,  . 

r,h  r,h 

Its  length  is  certainly  smaller  than  cr,  where  the 
constant  c   depends  only  on  the  smallest  angle  of 
inclination  of  a  tangent   y  to  the  x-axis. 

Between  the  values  of  v  at  P,   and  T,   we 
have  the  relation 


Squaring  both  sides  and  applying  Schwarz's  In- 
equality, we  obtain 

v(Pj^)^  <  2v(Pj^)^  +  2cr  •  h  V"  v^ 

Pr,h 
Summing  with  respect  to  P.   in  the  x-direction, 
we  get 


h  \  v^  <  2crv(T^)^  +  2c^r^h  V  v^    .      ' 


y  v^  <  2crv(Tj^)^  +  2c^r^h  V 
Pr  Pr 


Summing  again  in  the  y-direction  we  obtain  the 
relation 
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(iii)    ^  y  y"  ^^  -  ^°^  Y  ^^^h^^  "^  2c^r^h  y  V 


Writing  dovm  the  inequalities  associated  with 

the  other  portions  of   |    and  adding  all  the  inequal- 

19 

ities  together  we  obtain  the  desired  inequality   (13)  • 

We  next  set  'Vv,  =  ^v,  -  ^v  so  that  v,  =  0  on  \, 

Then  since  h'v~  y   (v  +  v  )   remains  boxmded 

h 
as  h  — >  0,  we  obtain  from  (13) 

2  ' 

(16)  ^^^v2<Xr 

r,h 

where  X  is  a  constant  which  does  not  depend  on  the 

function  V  or  the  mesh  size.   Extending  the  sum  on 

the  left  to  the  difference   S^.  ,  -  S^  .   of  two  bound- 

v,h    p,h 

ary  strips,  the  inequality   (16)   stills  holds  with 
the  same  constant  X  a-^d  we  can  pass  to  the  limit 
h  — >  0. 

From  the  inequality   (16)   we  then  get 
f 


1 

r 


2 
V  dxdy  <  )ic,  V  =  u-f 


S  -S^ 
r  p 


19 

By  similar  reasoning  we  can  also  establish  the  inequality 


(15)     h^  ^v^  <  c^h  y  v^  +  ^2^^  Y}""^ ""  ""P 

\  'k  \ 

in  which  the  constants  c,  and  c^  depend  only  on  the  re- 
gion G  and  not  on  the  mesh  division. 
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Now  letting  the  narrower  boundary  strip  S   approach 
the  boundary  we  obtain  the  Inequality 


^  I  j  v^dxdy  =  ^  I  (u-f  )^dxdy  <  Hi* 

S  S 

r  r 

which  states  that  the  limit  function  u  satisfies 

the  prescribed  boundary  condition. 

l^..   Applicability  of  the  method  to  other  problems 
Our  method  is  brsed  essentially  on  the  inequal- 
ities arising  from  the  lemmas  stated  previously  since 

the  principal  points  of  the  proofs  follow  from  the  two 

20 
last  lemmas  in  part  1  of  Section  i^»  No  use  Is  made 

of  special  fundamental  solutions  or  special  proper- 
ties of  the  difference  expression,  and  therefore  the 
method  can  be  carried  over  directly  to  the  case  of  ar- 
bitrarily many  independent  variables  as  well  as  to  the 
eigenvalue  problem, 

C^x^   ^y2 

The  same  sort  of  convergence  relations  will  obtain  in 

21 
this  case  as  above.    Also  the  method  applies  to 

20 

For  an  application  of  corresponding  Integral  inequal- 
ities see  K.  Priedrichs,  "Die  Rand-und  Eigen  Wert  prob- 
lems aus  der  Theorle  der  elastischen  Flatten",  Math.  Arnii 
98,  p.  222. 

21 

Similarly  one  proves  at  the  same  time  that  every  solu- 
tion of  such  a  differential  equation  problem  has  deriva- 
tives of  every  order. 
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linear  partial  differential  equations  of  other  types, 
in  particular  its  application  to  equations  with  varia- 
ble coefficients  requires  only  some  minor  modifications. 
The  essential  difference  in  this  case  lies  only  in 
proving  the  boundedness  of  h  ^  5Z  %  which  of  course 
does  not  always  hold  for  an  arbitrary  linear  problem. 
However  in  case  this  sum  Is  not  bounded  It  can  be 
shown  that  the  general  bo\mdary  value  problem  for  the 
differential  equation  in  question  also  possesses  ef- 
fectively no  solutions,  but  that  in  this  case  there 

exist  non-vanishing  solutions  of  the  corresponding 

22 

homogeneous  problem,  i.e.  elgenf unctions   • 

5.   The  boundary  value  problem  ZN  //^^  =  0 
In  order  to  show  that  the  method  can  be  carried 
over  to  the  case  of  differential  equations  of  higher 
order,  we  will  treat  briefly  the  boundary  value  problem 
of  the  differential  equation: 


<9x^     Jx2^y2   ^y^ 


We  seek.  In  G,   a  solution  of  this  equation  for 
which  the  values  of  u  and  Its  first  derivative  are 
given  on  the  boundary,  being  specified  there  by  some 

function  f(x,y). 

22 

See  Courant-Hllbert,  Methoden  der  mathematlschen  Physlk  1, 

Chap.  Ill,  Section  3»  where  the  theory  of  Integral  equations 
Is  handled  with  the  help  of  the  corresponding  principle  of 
the  alternative.   See  also  the  Dissertation  (Gottlngen)  of 
W.  V.  Koppenfels  which  will  appear  soon. 
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To  this  end  we  assume  as  previously  that  f(x,y) 
together  with  its  first  and  second  derivatives  is  con- 
tinuous in  that  region  of  the  plaxie  containing  the  re- 
gion G, 

We  replace  our  differential  equation  problem  by 
the  new  problem  of  solving  the  difference  equation 
/\u  =  0  in  the  mesh  region  G  subject  to  the  condi- 
tion that  at  the  points  of  the  boundary  strip  \]^'*'    I  ^ 
the  function  u  takes  on  the  values  f(x,y).   From 
Section  2  we  know  that  this  boundary  value  problem  has 
a  unique  solution.   We  will  show  that  as  the  mesh  size 
decreases  this  solution,  in  each  interior  subregion  of 
G,   converges  to  the  solution  of  the  differential  equa- 
tion, and  that  all  of  its  difference  quotients  converge 
to  the  corresponding  partial  derivatives. 

We  note  first  that  for  the  solution  u  =  u,  ,   the 

sum 

2t—  X-   ,..^     ^   0.-2   ^  2 


^TL  '"-  ^  '^U  *  v: 


remains  bounded  as  h  — >  0,   That  is,  by  applying  the 
mlnimvun  requirements  on  the  solution  (part  2,  Section 
2)  one  finds  that  this  s\ini  Is  not  larger  than  the  cor- 
responding sum 

h^  V  V  (f^  +  2f^   +  f^  ) 

2_tZ-   XX     xy    yy'    * 

^h 

and  this  converges  as  h  — >  0  to 
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^         2^     ,Z^  .2 


(cL|.  2^111.  +^)  dxdy 


which  exists,  by  hypothesis. 

From  the  boundedness  of  the  sxm. 


^^.Z  ^''^-  ■"  ^''"-  ■"  """-^ 


XX     xy    yy' 
follows  immediately  the  boundedness  of  h^  \~  \~(y\u)' 


°h 


2  T-  T— ,  »..^2 


■^  i;,!'^" 


and  also  that  of  » 

h^  y  V"  (u^  +  u^)   and  li^  y  ^  ^^   • 

^h  ^ 

That  is,  for  arbitrary  w  the  following  inequal- 
ity  holds  (see  footnote  19), 

(15)     h^^^w^  <  ch^  J^(w^  +  w^)  +  chyw^ 
^h  ^h  Th 

Then  if  one  substitutes  the  first  difference  quotients 
of  w  for  w  itself  in  this  inequality  and  applies 
the  expression  over  the  subregion  of  G,   for  which 
they  are  defined,  then  one  finds  the  further  inequality, 

where  again  the  constant  c  Is  independent  of  the  func- 
tion and  of  the  mesh  size.   We  apply  this  inequality  to 
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w  =  u,   and  thus  find  the  boundedness  of  the  siam  over 
h 

P,  +  P,   on  the  right-hand  side;  (by  definition  these 
boundary  sums  converge  to  the  corresponding  integral 
containing  f(x,y).   Thus  from  the  boundedness  of 


^h 


h^  Y"(u^  +  2u^  +  u^  ) 
2_,2_  XX     xy    yy' 

follows  the  boundedness  of  h^~  y(u^  +  u  )   and  thence 

that  of  h^  ^u^   . 

^h 
For  the  third  step  we  substitute  one  after  the 

other  the  expressions  /\u.  /\u  ,  /\u  ,  /\u   »  •••  , 

X     y     XX 

for  w  in  the  inequality 

G*"-  G^  G* 

(see  part  2,  Section  i|.)   where  G   is  a  subreglon  of 
G  containing  G    in  its  interior.   The  expressions 
introduced  all  satisfy  the  equation  /\w  =  0,   It  fol- 
lows then  that  for  each  expression  in  turn  and  for  all 

■ft  ^  p  p 

subregions      G       of     G     that  the    sums,      h  /    ^>    (w     +  w   ), 

that  is ,  h^^X^ Au^  +  Au^ ) ,  ti^Z^ j:^ Au^  +  Au^y ) , . . . 

are  bounded  together  with  the  sums: 

h^  J^u^,  h^^  ^(u^  +  u^),  and  h^  ^(A^i)^, 
^h        ^h  \ 

which  are  already  known  to  be  bounded. 

'  kk  -  - 
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Finally  we  substitute  into  the  Inequality   (10), 

in  place  of  w,   the  sequence  of  functions   u   ,u   ,u   , 

A^  A^   ^y 

u    ,••..,   for  which 


(Aw)^   ,   i.e.    hV  yCzXu   )2,  ... 


h  n 


are  botmded  as  shown  above.  We  then  find  that  for 
all  subregions  the  sums 

\  °h 

remain  bounded. 

Prom  these  facts  we  can  now  conclude  as  previously 
that  from  our  sequence  of  mesh  functions  a  subsequence 
can  be  chosen  which  in  each  Interior  subregion  of  G 
converges  (together  with  all  its  difference  quotients) 
\inlformly  to  a  limit  function   (  or  respectively  its 
derivatives)  which  is  continuous  in  the  interior  of  G. 

We  have  yet  to  show  that  this  limit  function  which 
obviously  satisfies  the  differential  equation  ^  /\u  =  0 
also  takes  on  the  prescribed  boundary  conditions.   For 
this  purpose  we  say  here  only  that,  analogous  to  the 
foregoing,  the  expressions 


(u-f )^dxdy  <  cr^  , 


Ifv2 


dxdy  <  cr' 


S 
r 


-  kS  - 
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hold^.   That  the  limit  function  fulfills  these  con- 
ditions may  be  seen  by  carrying  over  the  treatment  in 
part  3*  Section  4  to  the  function  u  and  its  first 
difference  quotient. 

From  the  uniqueness  of  our  boundary  value  problem 
we  see  furthermore  that  not  only  a  selected  subsequence, 
but  also  the  original  sequence  of  functions   u   pos- 
sesses the  asserted  convergence  properties. 


-^That  the  boundary  values  for  the  function  and  its 
derivatives  actually  are  assumed  is  not  difficult  to 
prove.   See  for  instance  the  corresponding  treatment 
of  K.  Friedrichs  (loc.  cit.). 


-  i+6  - 
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II.   The  Hyperbolic  Case 

Section  1,   The  Equation  of  the  Vibrating  String. 

In  the  second  part  of  this  paper  we  shall  consider 
the  Initial  value  problem  for  linear  hyperbolic  partial 
differential  equations.   We  shall  prove  that  vinder  certain 
hypotheses  the  solutions  of  the  difference  equations  con- 
verge to  the  solutions  of  the  differential  equations  as 
the  mesh  size  decreases. 

We  can  discuss  the  situation  most  easily  by  consider- 
ing the  example  of  the  approximation  to  the  solution  of 
the  wave  equation 

<if<i  .ifa  =  0 

We  limit  oxirselves  to  the  particular  initial  value  problem 
where  the  value  of  the  solution  u,   and  its  derivatives 
are  given  on  the  line   t  =  0,        , 

In  order  to  find  the  corresponding  difference  equation, 
we  construct  in  the   (x,t) -plane  a  sqxiare  grid  with  lines 
parallel  to  the  axes  and  with  mesh  width  h.   Following  the 
notation  of  the  first  part  of  the  paper  we  replace  the  dif- 
ferential equation  (1)  by  the  difference  equation 

^tt  -  ^x5E  =  °   • 
If  we  select  a  grid  point,   P  ,   then  the  corresponding  dif- 
ference equation  relates  the  value  of  the  function  at  this 
point  to  the  values  at  four  neighboring  points.   If  we  char- 
acterize the  four  neighboring  values  by  the  four  indices 
1,  2,  3,  k     (cf.  Fig.  5),  then  the  difference  equation 
assumes  the  simple  form 

-  hi  - 
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(2)  ^1  +  ^3  -  ^2  -  \  =  °   • 

Note  that  the  value  of  the  function  u  at  the  point   P 
does  not  appear  itself  In  the  equation. 

S 

• 

1 

+    .    + 

•      •       •      • 

3 

«    +    »+        •»••• 

+   .   +   .         .... 

Pig.   5.  Fig,   6. 

We  consider  the  grid  split  up  into  two  subgrids, 
indicated  in  Pig«  5  by  dots  and  crosses  respectively. 
The  difference  equation  connects  the  values  of  the  func- 
tion over  each  of  the  subgrids  separately,  and  so  we 
shall  consider  only  one  of  the  two  grids.  As  initial 
condition  the  values  of  the  function  are  prescribed  on 
the  two  rows  of  the  grid,   t  =  0  and  t  =  h.   We  can 
give  the  solution  of  this  initial  value  problem  expli- 
citly; that  iSj  we  express  the  value  of  the  solution 
at  any  point  S  in  terras  of  the  values  given  along 
the  two  initial  rowso   One  can  see  at  once  that  the 
value  at  a  point  of  the  row   t  =  2h  is  uniquely  de- 
termined by  only  the  three  values  ac  the  points  close 
to  it  in  the  two  first  rows^   The  value  at  a  point  of 
the  fourth  row  is  uniquely  determined  by  the  values 
of  the  solution  at  three  particular  points  in  the  sec- 
ond and  third  rows,  and  through  them  it  is  related  to 
certain  values  in  the  first  two  rowso   In  general  to  a 
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point  S  there  will  correspond  a  certain  region  of 
dependence  In  the  first  two  rows;  it  may  be  found  by 
drawing  the  lines  x  +  t  =  const,   and  x  -  t  =  const., 
through  the  point   S  and  extending  them  until  they 
meet  the  second  row  at  the  points   a  and  B  respec- 
tively (cf.  Pig.  6).   The  triangle   Sap  is  called  the 
triangle  of  determination  because  all  the  values  of  u 
in  it  remain  unchanged  provided  the  values  on  the  first 
two  rows  of  it  are  held  fixed.   The  sides  of  the  tri- 
angle are  called  lines  of  determination. 

If  one  denotes  the  differences  of  u  in  the  di- 
rection of  the  lines  of  determination  by  u  and  u, 
that  is, 

^1  =  ^1  -  %   '    ^1  '  ^1  "  ^2 

U2  =  U2  -  U3   ,    u^^  '  \'  ^2 

then  the   difference   equation  assvimes   the   form 

/  / 

u^  =  U2 

I.e.  along  a  line  of  determination  the  dif- 
ferences with  respect  to  the  other  direction  of  de- 
termination are  constant,  and  thus  are  equal  to  one  of 
the  given  differences  between  the  value  at  two  points 
on  the  first  two  rows.  Moreover  the  difference  u,,  -  u^ 
is  a  sum   of  differences  u   along  the  determining  line 
aS,  so  that  using  the  remark  above,  we  can  obtain  the 

-  hS   - 
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final  result  (in  obvious  notation): 

/ 


Pi 


(3)  ^S  =  ^a-'I]^ 

^1 

We  now  let  h  go  to  zero,  and  let  the  prescribed 
values  on  the  second  and  first  rows  converge  uniformly 

to  a  twice  continuously  differentiable  function,   f(x), 

/ 
and  the  difference  quotients  -^r^     there  converge  uni- 
formly to  a  continuously  differentiable  function  g(x). 
Evidently  the  right-hand  side  of   (3)   goes  over  uni- 
formly to  the  expression 


ik)  f(x-t)  +^   \   g(4)dC 

x-t 
if  S  converges  to  the  point   (x,t).   This  is  the  well- 
known  expression  for  the  solution  of  the  wave  equation 
(1)   with  initial  values  u(x,0)  =  f(x)   and  <^^i^'^^  = 

f'(x)  + /Z   g(xr).   Thus  it  is  shown  that  as  h >  0 

the  solution  of  the  difference  equation  converges  to 
the  solution  of  the  differential  equation  provided  the 
initial  values  converge  appropriately  (as  above). 

Section  2»   On  the  Influence  of  the  Choice  of  Mesh. 

The  Domains  of  Dependence  of  the  Difference  and  Differential 

Equations. 

Th,  relationships  developed  in  Section  1  suggest 

the  following  considerations. 

-  50  - 
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In  the  same  way  that  the  solution  of  a  linear  hy- 
perbolic equation  at  a  point  S  depends  only  on  a 
certain  part  of  the  initial  line  -  namely  the  "domain 
of  dependence"  lying  between  the  two  characteristics 
drawn  through  S,   the  solution  of  the  difference 
equation  has  also  at  the  point   S  a  certain  domain 
of  dependence  defined  by  the  lines  of  determination 
drawn  through  S,   In  Section  1  the  directions  for 
the  lines  of  determination  of  the  difference  equation 
coincided  with  the  characteristic  directions  for  the 
differential  equation  so  that  the  domains  of  depend- 
ence coincided  in  the  limit.   This  result  however 
was  based  essentially  on  the  orientation  of  the  mesh 
in  the   (x,t) -plane,  and  depended  furthermore  and  the 
fact  that  a  square  mesh  had  been  chosen.   We  shall 
now  consider  a  more  general  rectangular  mesh  with  mesh 
size  equal  to  h  (time  interval)   in  the  t-direotlon 
and  equal  to  kh  (space  interval)  in  the  x-direction, 
where  k  is  a  constant.   The  domain  of  dependence 
for  the  difference  equation,   u. r  -  u  -  =  0  for  this 
mesh  will  now  either  lie  entirely  within  the  domain 

of  dependence  of  the  differential  equation, 

2     2 
°— 3  -  «--5  =  0,   or  on  the  other  hand  will  contain  this 

latter  region  inside  its  own  domain  according  as 

k  <  1  or  k  >  1  respectively. 
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Prom  this  follows  a  remarkable  fact:   if  for  the 

case   k  <  1   one  lets   h  >  0,   then  the  solution 

to  the  difference  equation  in  general  cannot  converge 

to  the  solution  of  the  differential  equation.   In  this 

case  a  change  in  the  initial  values  of  the  solution  of 

the  differential  equation   in  the  neighborhood  of  the 

endpoints   a  and   (3   of  the  domain  of  dependence 

( see  Fig.  7. ) 

S 
causes  according  ^'^^ 


to  formula   (i|),  a  ^z'.   .   .  s 

'I   ....  \ 

change  in  the  solu-  /  \ 

/  \ 

tlon  itself  at  the        -^ ^^ 

point      (x,t).     For  / ••••       \ 

the  solution  of  the 


> 


difference  equation  „.    „ 

rig.  f. 

at  the  point  S  how- 
ever the  changes  at  the  points  a  and  (3  are  not  rel- 
evant since  these  points  lie  outside  of  the  domain  of 
dependence  of  the  difference  equations.  -  That  converg- 
ence does  occur  for  the  case  k  >  1  will  be  proved  in 
Section  3«   See  for  example  Fig,  9. 

If  we  consider  the  differential  equation 

(5)         2  jU  -  c^  -  ^  =  0 

dt^    dx^    dr 

in  two  space  variables,   x  and  y,   and  time,   t,   and 

-  SZ   - 
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if  we  replace  It  by  the  cori'e spending  difference  equa- 
tion on  a  rectilinear  grid,  then  in  contrast  to  the  case 
of  only  two  independent  variables  it  is  Impossible  to 
choose  the  mesh  division  so  that  the  domain  of  depend- 
ence of  the  difference  and  differential  equations  co- 
incide, since  the  domain  of  dependence  of  the  difference 
equation  is  a  quadrilateral  while  that  of  the  differ- 
ential equation  is  a  circle.   Later  (cf.  Section  L(.)  we 
shall  choose  the  mesh  division  so  that  the  domain  of 
determination  of  the  difference  equation  contains  that 
of  the  differential  equation  in  its  interior,  and  shall 
show  that  once  again  convergence  occurs. 

On  the  whole  an  essential  result  of  this  section 
will  be  that  in  the  case  of  each  linear  homogeneous  hy- 
perbolic equation  of  second  order  one  can  choose  the 
mesh  so  that  the  solution  of  the  difference  equation 
converges  to  the  solution  of  the  differential  equation 
as   h >  0,  (see  for  instance  Sections  3j  k->    7,  8). 

Section  3.   Limiting  Values  for  Arbitrary  Rectanp^ilar  Grids. 
Now  we  consider  further  the  wave  equation 

(1)  d!u  .cLu  =   0   , 

but  impose  it  now  on  a  rectangular  grid  with  time  in- 
terval h  and  space  interval  kh.  The  corresponding 
difference  equation  is 
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(6)   L(u)  =  ^in^-2u^^u^)    -  1   (u2-2u^+u^}  =  0 


where  the  indices  represent  a  "fundamental  rhombus" 

with  midpoint   P   and  corners  P,  ,  Pp*  P-j,  Pi   (see 

Fig.  8),  According 

S 


l^.    0    2 


Pig.  8 


•  •  • 


•  ;  •  • 


Fig.  9 


to  the  equation  L(u)  =  0  the  value  of  the  fxmction 
u  at  a  point  S  can  be  represented  by  its  values  on 
that  section  of  the  initial  rows  t  =  0  and  t  =  h 
cut  out  by  lines  of  determination  through  S  parallel 
to  the  sides  of  an  elementary  rhombus.  We  assume  that 
the  initial  values  are  prescribed  in  such  a  way  that 
as  t  — >  0  for  fixed  k  the  first  difference  quo- 
tients formed  from  them  converge  \inlformly  to  given  con- 
tinuous functions  on  the  line  t  =  0,  The  initial 

values  can  be  used  to  form  ajn.   explicit  representation 
of  the  solution-  of  the  difference  equation  (correspond- 
ing to  (3)  in  Section  1);  however  it  is  too  compli- 
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cated  to  yield  a  limiting  value  easily  as  h  — ->  0. 
Thus  we  will  try  another  approach  which  will  also  make 
it  possible  for  us  to  treat  the  general  problem,  ^ 

We  multiply  the  difference  expression  L(u)  by 
(u,-u,)   and  form  the  product  using  the  following 
identities: 


> 


(7)     {n^-\x^){vL^-2^i^+n^)   =  (u^-u^)^  -  (u^-u^)^ 


(8)     (u^-U3)(u2-2u^+u^)  =  (uj^-Uq)^  -  (u^-u^)^ 


-  ■|[{u^-U2)^  +  (u^-u,  )^  -  (u2-u^)^  -  (u, -u^)^]   . 


Then  we  obtain 


(9)      2(u^-U3)L(u)   =^2^1--^)    [^^1"%)^  -   ^%-^3)^l   + 

*7r;7    [^^1-^2^^  *   ^^1-^4)^  -   (U2-U3)2  -   (u^-U3)2] 

The  product  (9)  Is  now  summed  over  all  elementary 
rhombuses  of  the  domain  of  determination.   Sap,  The 
quadratic  difference  terms  on  the  right-hand  side  al- 
ways appear  with  alternate  signs  in  two  neighboring 
rhombuses  so  that  they  cancel  out  for  any  two  rhom- 
buses belonging  to  the  triangle   Sap,   Only  the  s\ams 


^For  the  following  see  K,  Priedrichs  and  H.  Lewy,  "Uber 
die  Eindentigk:eit..,etc.,  Math.  Ann.  (98,  1928)  p.  192  ff, 
where  a  similar  transformation  is  used  for  integrals. 
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of  squared  differences  over  the  "boundary"  of  the  tri- 
angle remain,  and  we  obtain  the  relation: 

Sap  Sa     ^ 

(10)  +  hy[2(l-^)(|)^  +  ^2^^'  )^]  + 

I  II      ^  k  1£ 

Here  u   and  u  denote  differences  in  the  direction 
of  determination  defined  in  Section  1,  while  t1  des- 
ignates the  difference  of  the  functional  values  at  two 
neighboring  points  on  a  mesh  line  parallel  to  the  t 

axis.   The  range  in  rr-     over  which  (u')    is  taken 

Sa 
is  the  outermost  boundary  line   Sa  and  its  nearest 

parallel  neighbor  fo\ind  by  shifting  the  points  of  Sa 

downward  by  the  amoxuit  h.   There  is  a  similar  range 

for   (u  )    in  YZ*      ^"^  3°  *^^  °^  ^^®  differences, 

sp 

u^,  u  ,   and  il  appear  once  and  only  once. 

For  the  solution  to  the  problem  L(u)  =  0  the 
.right-hand  side  of  (10)   vanishes.   The  sum  over  the 
initial  rows   I  and  II  which  occurs  there  remains 
bounded  as  h  — >  0  (for  fixed  k)j   specifically  it 
goes  over  into  an  integral  of  the  prescribed  function 
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along  the  initial  line.   Consequently  the  sums  over 

Sa  and   Sp  in   (10)   also  remain  bounded.   If  now 

k  >  1  as  we  must  require  (see  previous  discussion), 

then  1  -  — »j  is  non  negative,  and  we  find  that  the 

k 
individual  sums 

Sa  Sp 

extendt.  d  over  any  line  of  determination  whatever,  re- 
main bounded. 

Prom  this  we  can  derive  the  "uniform  continuity" 
(equi-continuity)  (cf.  Section  I4.   of  the  first  part  of 
the  paper)  of  the  sequence  of  grid  functions  in  all 
directions  in  the  plane j    since  the  values  of  u 
on  the  initial  line  are  bounded,  there  must  exist  a 
subset  which  converges  uniformly  to  a  limit  function 
u(x,t ) . 

Both  the  first  and  second  difference  quotients 
of  the  function  u  also  satisfy  the  difference  equa- 
tion L(u)  =  0.   Their  initial  values  can  be  expressed 
through  the  equation  L(u)  =  0  in  terras  of  the  first, 

second  and  third  difference  quotients  of  u  Involving 

— . — ^ — >_ 

-^If  S-  and  Sp  are  two  points  at  a  distance   6,   then 
one  connects  tnem  by  a  path  of  two  segments,  Sj^S  and  SSg, 
vAiere  the  former  is  parallel  to  one  line  of  determination 
and  the  latter  to  the  other.   Then  one  finds  the  appraisal. 


)2 


S^S         ^   SS2 


-  57  - 


NYO-7689 


initial  values  at  points  on  the  two  initial  lines   I 
and  II  only.   We  require  that  they  tend  to  continu- 
ous limit  functions,  that  is  that  the  given  initial 
values  u(x,0),  u.(x,0)  be  three  times  or  respec- 
tively twice  continuously  differentiable  with  respect 
to   X. 

From  this  we  can  apply  the  convergence  consider- 
ations set  forth  above  to  the  first  and  second  dif- 
ference quotients  of  u,   as  well  as  to  u   Itself, 
and  we  can  choose  a  subsequence  such  that  these  dif- 
ference quotients  converge  uniformly  to  certain  func- 
tions, which  must  be  the  first  or  respectively  second 
derivatives  of  the  limit  function  u(x,t).   Hence  the 

limit  function  satisfies  the  differential  equation 

2     2 
<i— ^  -  <i-^  =  0  which  results  as  the  limit  of  the  dif- 

ference  equation  L(u)  =  0;   it  represents  indeed  the 

solution  of  the  initial  value  problem.   Since  such  a 

solution  is  viniquely  determined,  every  subsequence  of 

mesh  fvmctions  converges,  and  therefore  the  sequence 

itself  converges  to  the  limit  function. 

Section  li.   The  Wave  Equation  in  Three  Variables 
We  treat  next  the  wave  equation, 

(11)       2fltfu   A  .cifu^  0 
dt^  ^x^       <^7 

and  consider  its  relation  to  the  observations  on  the 
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domain  of  dependence  discussed  in  Section  2,   The  do- 
main of  dependence  of  the  differential  equation   (11) 
is  a  circular  cone  with  axis  parallel  to  the  t-direc- 
tion  and  with  apex  angle   a,   where   tan  a  =/2«    "^^ 
any  rectilinear  grid  parallel  to  the  axes  we  introduce 
the  corresponding  difference  equation 


^12)         2-tl   -  ^x5  -  V  =   ° 


This  equation  relates  the  functional  values  of  u  at 
points  of  an  elementary  tetrahedron.   In  fact  it  allows 
the  value  of  the  function  u  at  a  point  S  to  be  ex- 
pressed uniquely  in  terms  of  the  values  of  the  function 
at  certain  points  of  the  1"-/o  initial  planes  t  =  0  'and 
t  =  h.  At  each  point  S  we  obtain  a  pyramid  of  deter- 
mination which  cuts  out  from  the  two  oase  planes  two 
rhombuses  as  domains  of  dependence. 

If  we  let  the  mesh  widths  tend  to  zero,  keeping 
their  ratios  fixed,  then  we  can  expect  convergence  of 
the  sequence  of  mesh  functions  to  the  solution  of  the 
differential  equation  only  provided  the  pyramid  of 
determxnation  contains  the  cone  of  determination  of 
the  differential  equation  in  its  interior.   The  simp- 
lest grid  with  this  property  will  be  one  constructed 
in  such  a  way  that  the  pyramid  of  determination  Is 
tangent  to  the  exterior  of  the  cone  of  determination. 
Our  differential  equation  is  chosen  so  that  this  occurs 
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for  a  grid  of  cubes  parallel  to  the  axes. 

The  difference  equation  (12),   In  the  notation 
of  Fig.  10,   assumes  for  such  a  grid  the  form: 

(13)   L(u)  =  72^^a'-2^o-^a^  "  h^^'l'^^'o'^^^    "" 
h  n 

n.        ^ 

in  which  the  functional  value,  u  ,  at  the  midpoint 
actually  cancels  out.   The  values  of  the  solution  on 
the  two  Initial  planes  must  be  the  values  of  a  function 
having  fotir  continuous  derivatives  with  respect  to 

^1'  ^2*^V 

For  the   convergence   proof 

we  again  use  the  method  de- 
veloped in  Section  3.      We 
construct  the   triple   sum 


u. 


.2t—  t-  t-o   "a 


/. 


FigT  10. 


XEP 


u. 


-S  L(u)   =   0 


for  the  solution  to  the  difference  equation,  where  the 
summation  is  to  be  extended  over  all  elementary  octahe- 
drons of  the  pyramid  of  determination  emanating  from 
the  point  S,  Then  almost  exactly  as  before  we  find 
that  the  values  of  the  function  u  at  the  Interior 
points  of  the  pyramid  of  determination  cancel  out  in 
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the  aimmatlon  and  that  only  the  values  on  the  two 
pyrajnlds  called  F,  and  on  the  two  base  surfaces 
I  and  II  remain. 

If  we  denote  by  u'   the  difference  of  the  values 
of  the  function  at  two  points  connected  by  a  line  of 
an  elementary  octahedron,  then  we  can  write  the  result 
as 

ilk)   X!II  ^""'^^  "HE  ^""^^^  "  °    ' 

F  I   II 

where  the  sum  Is  extended  over  all  surfaces  containing 

differences  u';  each  such  difference  is  to  appear 

26 
only  once.    The  double  sum  over  the  two  initial  sxir- 

faces  stays  bounded  since  it  goes  over  into  an  inte- 
gral of  the  initial  values.   Therefore  the  sum  over 
the  'burface  of  determination"  F  remains  bounded. 
We  now  apply  these  results  not  to  u  itself, 
but  to  its  first,  second  and  third  difference  quotients, 
which  themselves  satisfy  the  difference  equation  (13 )• 
Their  initial  values  can  be  expressed  using  only  values 
on  the  first  two  initial  planes  by  means  of   (13)  using 

first  through  fourth  difference  quotients.   If  w  =  w^^ 

■Tg — ' — ~ — 

The  grid  ratio  has  been  chosen  in  such  a  way  that  the 
differences  between  values  of  u  appearing  on  the  two 
neighboring  surfaces  in  F  do  not  occur,  (as  they  did  in 
the  general  case  in  one  dimension  treated  in  Section  3)« 
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is  one  of  the  difference  quotients  of  any  order  up  to 
third  order,  then  we  know  that  the  sum  h  >  T"(^) 
extended  over  a  surface  of  determination  remains  bovmd- 
ed.   Prom  this  it  follows,  through  exactly  the  methods 
used  in  Section  ij.  of  the  first  part  of  the  paper,  that 
the  function  u  and  its  first  and  second  difference 
quotients  are  uniformly  continuous  (equi-continuous) . 
Thus  there  exists  a  sequence  of  mesh  widths  decreasing 
to  zero  such  that  these  quantities,  which  are  bounded 
initially,  converge  to  continuous  limit  functions,  and 
in  fact,  converge  to  the  solution  of  the  differential 
equation  and  to  the  first  and  second  derivatives  of 
this  solution,  all  exactly  as  we  found  earlier  (Section  3)» 
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APPENDIX 
Supplements  and  Generalizations 

Feetion  5«   Example  of  a  Differential  Equation  of  First  Order 

We  have  seen  in  Section  2  that  in  the  case  when 
the  region  of  dependence  of  the  difference  equation 
covers  only  a  part  of  the  region  of  dependence  of  the 
differential  equation,  the  influence  of  the  rest  of  the 
region  is  not  included  in  the  liTn?'t.  We  can  demon- 
strate this  phenomenon  explicitly  by  the  example  of 
the  differential  equation  of  first  order. 

If  we  replace  it  by  the  difference  equation 
(15)  ^^t  "  ^x  *  ^x  *  °   • 


In  the  notation  of  Fig.  5  this  becornes 
u„  +  u, 


(16)         u,  = 


_  2 


1      2 

As  before,  the  difference  equation  connects  only  the 
points  of  a  submesh  with  one  another.   The  initial 
value  problem  consists  of  assigning  as  initial  values 
for  u  at  points  x  =  2ih  on  the  row  t  =  0  the 
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values,  fpi^'      assumed  there  by  a  continuous  function 
f(x). 

We  consider  the  point  S  at  a  distance   2nh  up 
along  the  t  axis.   It  is  easy  to  verify  that  the 
solution  u  at   S  is  represented  as 


(17)         U3=|^  TS;("<'  ' 


P^.x  n+i'  "2i 
i=-n  '^ 


As  the  mesh  size  decreases,  that  is  as  n  — >  oo  the 

sum  on  the  right-hand  side  tends  simply  to  the  value 

f  .   This  can  be  demonstrated  from  the  continuity  of 
o 

f(x)   and  from  the  behavior  of  the  binomial  coefficients 
as  n  increases!  (see  the  following  paragraph). 

Section  6.   The  Equation  of  Heat  Conduction. 

The  difference  equation  (16)  of  Section  5  can 
also  be  interpreted  as  the  analogue  of  an  entirely  dif- 
ferent differential  equation,  namely  the  equation  of 
heat  conduction, 

(le)  24t  -  fi^   =    0 

In  any  rectangular  mesh  with  mesh  spacing  ^     and 
h  in  the  time  and  space  directions  respectively  the 
corresponding  difference  equation  becomes 

u,  -  u        Uo  +  U|  -  2u^ 

(19)        2(-^--^)   =   (-^ 1 ^) 

J  h^ 
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In  the  limit  as  the  mesh  size  goes  to  zero  the  dif- 
ference equation  preserves  its  form  only  if  U     and 

2 
h   are  decreased  proportionately.   In  particular  if 

/     2 
we  set  ^  =  h  ,   then  the  value  u   drops  out  of 

the  equation  and  the  difference  equation  becomes 

Up  +  u, 
(16)  u,  =  -^ ^ 

^  Z 

The  solution  to   (16)   is  given  by  formula 


(17)       -(0.t)=f:-fe(n;i)f2i    • 

i»  -n 

As  the  mesh  width  decreases,   a  point     4     on  the 
X-axis   is  always  characterized  by  the    index 

(20)  2i  =  I 

The  mesh  width  h  is  related  to  the  ordinate  t  of 
a  particular  point  by 


(21)  2nh^  =  t 


We  shall  investigate  what  happens  to  formula  (17) 
as  h  — >  0,  that  is  n  — >  oa  Using  (21)  we  write 
(17)   in  the  form 


(22)       -(Q>^)=I  ,^^   (n;i)^2i'  2h 
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For  the  coefficient  of  Shfp^  =  2hf(i)   we  use 
the  abbreviation 

The  limiting  value  of  the  coefficient  which  is  usually 
determined  by  using  Stirling's  formula,  we  will  calcu- 
late here  by  considering  the  function  BprS^^      ^^   ^^® 
solution  of  an  ordinary  difference  equation  which 
approaches  a  differential  equation  as  h  -->  0,   As 
the  difference  equation  one  finds 

k  I gh(«*2h)  -  g^(c) )  =  -  k  Sh'^'  IWtt 

(in  which  we  have  written  gv(4)   instead  of  62n^^^^*   °^ 

■'I 

gv^(^)   satisfies  the  normalization  condition 


t+h4+2h^ 


^  e^i^)   •  2h  =  2/t 
i=-n 

This  siam  is  over  the  region  of  dependence  of  the  dif- 
ference equation,  and  as  h  — >  0  this  covers  the  en- 
tire X-axis. 
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It  can  be  shown  easily  that  gv>(4)   converges 
uniformly  to  the  solution  g(x)   of  the  differential 
equation 

g«(x)   =  -g(x)  I 

with  the  auxiliary  condition  '^ 

oo 

g(x)dx  =  2/t 

-00 

From  formula  (22)      after  passing  to  the  limit  we  find 


00 


2t 
^.,^     e  f(4)d4   , 

-00 


which  is  the  known  solution  of  the  heat  conduction 
equation. 

The  results  of  this  section  can  be  carried  over 
directly  to  the  case  nf  the  differential  equation. 


cit  d^        67 


and  so  on  for  even  more  independent  variables. 

Section  7«   The  General  Homoff;eneous  Linear  Equation  of 
Second  Order  in  the  Plane. 
We  consider  the  differential  equation 
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(23)  <^-^^^^a<i^  +  ^^^r^  =  o 

The  coefficients  are  assumed  to  be  twice  continuously- 
different  lable  with  respect  to  x  and  t,   while  the 
Initial  values  on  the  line  t  =  0  are  assumed  three 
times  continuously  differentlable  with  respect  to  x. 
We  replace  the  differential  equation  by  the  difference 
equation 

(24)  L(u)  =  u^^(x,t)-k^u^(x,t)+au^+Pu^+YU  =  0 

in  a  grid  with  time  mesh  width  h  and  space  mesh 
width  ch  so  that  in  a  neighborhood  of  the  appropri- 
ate part  of  the  initial  value  line  the  inequality 

2 
1  -  ^  >  e  >  0  holds  for  the  constant,   c.   The  ini- 

c 
tial  values  are  to  be  chosen  as  in  Section  3* 

For  the  proof  of  convergence  we  again  transform 
the  sum, 

P  U-,  -  U-j 

Sap        " 

by  using  identities   (7)   and   (8).   In  addition  to  a 
sum  (see  for  example  (10))  over  the  doubled  boundary 
of  the  triangle   Sap   (Fig.  6)   one  obtains  a  sum 
over  the  entire  triangle   Sap  because  of  the  varia- 
bility of  the  coefficient  k  and  the  presence  of 
lower  order  derivatives  in  the  differential  equation. 
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By  using  the  differentiability  of  k  and  the  Schwarz 
Inequality  one  can  show  that  this  latter  sum  is  bound- 
ed from  above  by 

Sap 

where  the  constant  C   is  independent  of  the  function 
u,   the  mesh  width   h,  and,  in  a  certain  neighborhood 
of  the  initial  line,  also  independent  of  the  point  S» 
Again  we  can  estimate  an  upper  bound  for 

Sap 

Sap         I  II 

where  the  same  properties  hold  for  C,   and  Cp  as 
are  stated  above  for  C, 

Thus  we  obtain  an  inequality  of  the  form 

hV  [2(1-4)(|)2  +i£^  (B')2  ] 


(25)     -H  hV  [2(1  -^)(§)2  +^  (^')^  ] 

sp      <*       e 

Sap  ^ 

■^ . _ 

'For  proof  one  makes  use  of  the  inequality  used  in 
footnote  25. 
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where   D,   for  all  points   S  and  mesh  widths  h, 
is  a  fixed  bound  for  the  sums  over  the  initial  line. 

As  vertices  of  our  triangles  we  choose  a 
sequence  of  points   S, ,  S, ,  ...  ,  S  =  S  lying  on  a 
line  parallel  to  the   t-axis.   By  summing  the  cor- 
responding sequence  of  inequalii;ies   (25)   after 
multiplying  by  h  we  obtain  the  following  inequality 

S  S  a        c         c   ^ 
o 


S  S^S        c         c 

0 

,  ,2     /  2     X  2 

<  nh^C^i:  H  [  (f )   +  (§  )   +  (g  )   ]  +  nhD 


Now  if  we  notice  that  one  can  express  a  dlf- 
ference  u'  or  u   in  terms  of  two  differences  t1 
and  a  difference  u   or  respectively  u',   then  we 
see  that  the  left-hand  side  of   (26)   is  larger  than 
the  simpler  form 

with  a  suitable  constant  C,  . 

Then  by  confining  the  discussion  to  a  suf- 
ficiently small  neighborhood,   0  <  t  <  nh  =  6   of 
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the  initial  line  where  5  is  small  enough  so  that 


C.  -  nhC^  =  C^  >  0 


we  find  from  ( 26 ) , 


o        '  ,^  2     /  2    „\  2     C, 

^   Sap     ^      ^       ^^        ^3 


The  bound  given  by   (2?)   when  combined  with 
(25)   gives  a  a  bound  on 


Sa  '^      Sp  " 


from  which,  as  in  Section  3»   one  can  prove  the 
viniform  continuity  of  u. 

We  apply  the  inequality  (25)  now,  instead 
of  to  the  function  u  itself,  to  its  first  and 
second  difference  quotients,  w,  which  also  satis- 
fy difference  equations  whose  second  order  terms 
are  the  same  as  those  of  ^2\\.)  *      In  the  rest  of  the 
terms  there  will  appear  lower  order  differences  of 
u  which  cannot  be  expressed  in  terms  w,  but  they 

will  appear  in  the  above  argximent  in  the  form  of 

2 

quadratic  double  sums  multiplied  by  h  .  This  is 

enough  to  let  us  reach  the  same  conclusions  for  the 
difference  equation  for  w  as  we  found  above  for 
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u.   So  we  can  conclude  from  this  the  uniform  contin- 
uity (equi-continuity )  and  boundedness  of  the  func- 
tion u  and  its  first  and  second  derivatives. 
Consequently  a  subsequence  exists  which  converges 
uniformly  to  the  solution  of  the  initial  value 
problem  for  the  differential  equation.  Again  from 
the  uniqueness  of  the  solution  we  find  that  the  se- 
quence of  functions  itself  converges. 

In  all  of  this  the  assiimption  must  be  made 
that  the  difference  quotients  up  to  third  order 

involving  values  on  the  two  initial  lines  converge 

2fi 

to  continuous  limit  functions. 

Section  8.   The  Initial  Value  Problem  for  an 
Arbitrary  Linear  Hyperbolic  Differential  Equation 

of  Second  Order. 
We  shall  now  show  that  the  methods  developed 
so  far  are  adequate  for  solving  the  initial  value 
problem  for  an  arbitrary  homogeneous  linear  hyper- 
bolic differential  equation  of  second  order.   It 
suffices  to  limit  the  discussion  to  the  case  of 
three  variables.   The  development  can  be  extended 
immediately  to  the  case  of  more  variables.   It  is 

This  assumption  and  also  the  assumptions  on  the 
differentiability  of  the  coefficients  of  the  dif- 
ferential equation,  and  further  on  the  restriction 
to  a  sufficiently  small  region  of  the  initial  line 
can  be  weakened  in  special  cases. 
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easy  to  see  that  a  transformation  of  variables  can  reduce 
the  most  general  problem  of  this  type  to  the  following: 
a  function  u(x,y,t)   is  to  be  found  which  satisfies  the 
differential  equation 

(28)     u^^  -  (au^  +  2bu^y  +  cUyy)  +  au^  +  pu^  +  yUy  +  6u  =  0, 

and  which,  together  with  its  first  derivative,  assumes  pre- 
scribed values  on  the  surface  t  =  0.   The  coefficients  in 
equation  (28)   are  functions  of  the  variables  x,y,   and 
t   and  satisfy  the  condition 

2 

a>0   ,   c>0   ,   ac-b  >0 

We  assume  that  the  coefficients  are  three  times  dlf- 
ferentiable  with  respect  to  x,  y,   end  t,   and  that  the 
initial  values  u  and  u.   are  respectively  four  and 
three  times  continuously  differentiable  with  respect  to 
X  and  7* 

The  coordinates  x  and  y  can  be  drawn  from  a  given 
point  on  the  initial  plane  in  such  a  way  that  b  =  0 
there.  Then  of  course  in  a  certain  neighborhood  of  this 
point  the  conditions 

a  -  |b|  >  0  ,   c  -  |b|  >  0 

hold.  We  restrict  our  investigation  to  this  neighborhood. 
Wc  can  choose  a  three  times  continuously  differentiable 
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function  d  >  0  so  that 


(29) 


for  some  constant  e.   Then  we  put  the  differential  equa- 
tion Into  the  form 


(30) 


u^o.  -  (a-d)u   -  (c-d)u   -  ^(d+b)(u   +  2u   +  u   ) 
"tt   ^   ^'  XX   ^^         77   ^  '^  XX   '^  xy    yy^ 

-  |(d-b)(u^  -  2u^  +  Uyy)  +  au^  +  pu^  +  yu^  +  5u  =  0 


We  now  construct  in  the  space  a  grid  of  points, 

t  =/h,  X  +  7  =  mkh,  x  -  7  =  nk..^  (^,ra,n  =  .,.-1,0,1,2...) 

and  we  replace  equation  (30)  b7  a  difference  equation 
L(u)  =  0  over  this  mesh.   We  do  this  b7  assigning  to  each 
point  P   the  following  neighborhood:  The  point  P  /  or 
the  point  P   which  lies  at  a  distance  h  or  -h  re- 
spectivel7  along  the  t-axis  from  P  ;  also  the  points 
P,,  .«.,  Pg  which  lie  in  the  some   (x,y) -plane  with  P 
(see  Pig.  11).   These  points  constitute  an  "eleinentar7  oct- 


ahedron" with  vertices   P^/ , 

a 

^a'  ^1'  ^2'  ^3*  ^i^.*  ^°^   ®^°^ 

grid  point  lying  in  the  interior  o 

of  G  we  replace  the  derlva- 

k 

tives  appearing  in  equation  (30)       * 

7 

by  difference  quotients  over  ', 


f 


->  X 


Fig.   11 


-  74  - 


NYO-7689 


the  elementary  octahedron  about   P  . 
We  replace 

^tt  ^y  h  (^a.-  2%  +  uj 

h 
u^^  by   g  o  (^o  -  2u  +  u,  ) 


u__  -  2u„_  -H  u    by  -^  (u^  -  2u^  +  u^) 


XX     xy     yy 


k"h' 


The  first  derivatives  in   (30)   are  replaced  by  the 
corresponding  difference  quotients  in  the  elementary  oct- 
ahedron.  The  coefficients  of  the  difference  equation  are 
given  the  values  assumed  by  the  coefficients  of  the  dif- 
ferential equation  at  the  point   P  . 

On  the  initial  planes   t  =  0  and  t  =  h  we  assxame 
that  the  values  of  the  function  are  assigned  in  such  a 
way  that  as  the  mesh  size  approaches  zero  for  fixed  k, 
the  function  approaches  the  prescribed  initial  value,  and 
the  difference  quotients  over  the  two  planes  up  through 
differences  of  fourth  order  converge  uniformly  to  the  pre- 
scribed derivatives. 

The  solution  of  the  difference  equation  Liu)  =  0 
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at  a  point  Is  uniquely  determined  by  the  values  on  the 
two  base  surfaces  of  the  pyramid  of  determination  pas- 
sing through  the  point. 

To  prove  convergence  we  construct  a  sum 


•^'IIZ^ 


a'   a 


L(u) 


over  all  the  elementary  octahedrons  of  the  pyramid  of 
determination,  and  we  transform  it  by  using  identities 
(7)   and   (8).   This  gives  one  space  summation  multi- 
plied  by  h   and  qxiadratic  in  the  first  difference 

quotients,  and  also  over  a  double  surface  a  sum  which  is 

2 

multiplied  by  h   and  contains  the  squares  of  all  the 

difference  quotients  of  the  type  u„-  u  ,  u^-  u, ,  •••«. 

"^aoai 

..,u  -  Uq.   In  this  expression  according  to   (29)   the 
coefficients  are  larger  than  some  fixed  positive  constant 
in  all  those  cases  where  the  ratio   /k  between  the 
time  and  space  mesh  sizes  is  taken  sufficiently  small. 

Prom  here  on  one  can  proceed  exactly  as  before 
(Sections  7,  i^.)   and  prove  that  the  solution  of  the  dif- 
ference equation  converges  to  the  solution  of  the  dif- 
ferential equation. 
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